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SUMMARY 


Sea  water  oscillating  across  the  Earth's  magnetic  field  induces 
eddy  currents  in  the  sea  that  generate  an  electromagnetic  field  charac¬ 
teristic  of  the  velocity  field  in  the  sea.  The  electromagnetic  field  is 
the  sum  of  an  electrostatic  type  field,  resulting  from  accumulation 
of  surface  charge,  and  a  transverse  electric  type  field  having  an  impe¬ 
dance  proportional  to  phase  speeds  of  ocean  waves,  which  are  no  greater 
than  a  few  hundred  meters/sec  for  surface  waves  and  of  the  order  of 
1  m/sec  or  less  for  internal  waves.  The  electrostatic  part  of  the  field 
vanishes  in  polar  regions  and  is  strongest  in  equatorial  regions.  Con¬ 
versely,  the  transverse  electric  part  of  the  field  is  strongest  in  polar 
regions  and  weakest  in  equatorial  regions. 

Magnitude  of  the  electrostatic  part  of  the  field  produced  by  a 
one  meter  high  surface  wave  progressing  across  magnetic  meridians 
in  equatorial  regions  increases  with  wave  frequency  for  wavelengths 
less  than  an  ocean  depth  and  reaches  a  lower  limiting  value  of  the  order 
of  a  few  microvolts/m  at  wavelengths  greater  than  an  ocean  depth. 
Above  the  surface,  the  electric  field  vector  is  circularly  polarized  in 
a  vertical  plane  containing  the  direction  of  wave  propagation  and  rotates 
in  an  opposite  sense  to  the  velocity  vector  of  the  wave.  Immediately 
below  the  surface,  polarization  is  elliptical  with  the  same  sense  of 
rotation  as  the  velocity  vector.  At  the  ocean  floor,  the  field  is  linearly 
polarized  in  a  vertical  direction,  so  that  the  electrostatic  part  of  the 
field  vanishes  below  the  ocean  floor. 
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Similarly,  the  electrostatic  part  of  the  field  produced  by  an 
internal  wave  vanishes  both  above  and  below  the  ocean,  as  the  surface 
displacement  resulting  from  an  internal  wave  is  negligible.  Polariza¬ 
tion  of  the  field  within  the  sea  is  elliptical  with  the  same  sense  of  rota¬ 
tion  as  the  velocity  vector  of  an  internal  wave.  Magnitude  of  the  field 
per  unit  wave  height  at  a  thermocline  increases  with  frequency  for 
wavelengths  less  than  a  thermocline  depth  and  reaches  a  lower  limit¬ 
ing  value  of  the  order  of  a  few  tenths  of  microvolts/m  “  at  wavelengths 
greater  than  a  thermocline  depth. 

Surface  magnitudes  of  both  the  electric  and  magnetic  compo¬ 
nents  of  the  transverse  electric  part  of  the  field  decrease  with  wave 
number  for  wavelengths  of  surface  and/or  internal  waves  less  than  an 
ocean  depth.  Field  strengths  per  unit  wave  height  produced  by  surface 
waves  in  polar  regions  attain  maximum  values  of  the  order  of  a  few 
microvolts/m  and  a  few  tens  of  gamma/m  for  wavelengths  much  greater 
than  an  ocean  depth.  Internal  waves  produce  maximum  field  strengths 
at  the  surface  in  polar  regions  of  the  order  of  10  microvolts/m  and 
0.1  gamma/m  at  wavelengths  of  the  order  of  an  ocean  depth.  In  equa¬ 
torial  regions,  maximum  field  strengths  are  less  than  polar  values  by 
about  a  factor  of  two,  for  waves  progressing  along  magnetic  meridians. 
The  transverse  electric  part  of  the  field  vanishes  for  waves  propagat¬ 
ing  across  magnetic  meridians  in  equatorial  regions. 

The  magnetic  vector  of  the  transverse  electric  part  of  the 
field  is  polarized  in  a  vertical  plane  containing  the  direction  of  wave 
propagation,  and  the  electric  field,  by  definition,  is  linearly  polarized 
in  a  horizontal  direction  transverse  to  the  direction  of  wave  propaga¬ 
tion.  Polarization  of  the  magnetic  field  below  the  surface  exhibits  a 
complicated  pattern  with  both  depth  and  latitude.  At  and  above  the  sur- 
iace,  however,  the  magnetic  field  is  always  circularly  polarized  with  a 
sense  of  rotation  opposite  to  that  of  the  velocity  field. 
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I.  INTRODUCTION 


As  Faraday  [1832]  first  recognized,  sea  water  moving  across 
the  Earth’s  magnlfcic  field  induces  an  electric  current  in  a  sea.  First 
investigations  of  the  effect  were  concerned  with  measurement  of  elec¬ 
tric  fields  induced  by  steady  ocean  currents.  Efforts  culminated  in 
development  of  a  so-called  geomagnetic  electrokinetographic  or  GEK 
method  [von  Arx,  1950]  of  inferring  velocities  of  bioad,  steady  ocean 
currents  from  measurements*  of  horizontal  gradients  in  the  electric 
potential  at  the  ocean  surface.  Longeut-Higgins,  et  al.  [1954]  develop 
a  sound  basis  for  interpreting  GEK  measurements. 

Although  Longeut-Higgins,  et  al.  [1954]  consider  electric  fields 
and  currents  produced  by  ocean  waves  as  well  as  steady  currents, 
investigation  of  magnetic  effects  of  either  steady  currents  or  ocean 
waves  has  been  largely  neglected  until  recently.  Crews  and  Futterman 
[1962]  are  among  the  first  to  investigate  magnetic  fields  produced  by 
ocean  waves.  They  neglect  magnetic  induction  effects,  however,  and, 
by  using  the  integral  form  of  Ampere's  law,  determine  magnetic  fields 
produced  above  the  ocean  surface  by  slowly  varying  currents  in  a  deep 
ocean.  Warburton  and  Caminiti  [1964]  extend  the  analysis  to  determine 
magnetic  fields  produced  below  the  surface  of  a  deep  ocean,  and 
Groskaya,  et  al.  [1972],  to  shallow  seas. 


*  Potential  gradients  are  measured  by  towing  electrodes  separated  by 
distances  of  a  few  tens  of  meters  behind  a  ship. 


A  quasistatic  approximation  (neglecting  magnetically  induced 
electric  fields)  is  perhaps  most  appropriate  for  determining  magnetic 
fields  induced  by  tidal  motions  where  vorticity  produced  by  Coriolis 
forces  is  important.  For  motions  having  periods  appreciably  less  than 
tidal  periods,  however,  flow  fields  are  effectively  irrotational,  and 
magnetic  induction  effects  are  readily  included. 

Weaver  [1965]  is  the  first  to  have  accounted  for  magnetic 
induction  effects  in  determining  magnetic  fields  induced  by  surface 
waves  on  a  deep  ocean.  Beal  and  Weaver  [1970]  also  investigate  mag¬ 
netic  fields  produced  by  internal  waves  propagating  along  a  sharp 
thermocline.  Both  papers,  however,  neglect  the  electric  field. 

Measurements  by  Maclure,  et  al.  [1964]  of  magnetic  fields 
induced  at  the  surface  by  ocean  swell  agree  with  estimates  made  by 
Weaver.  Measurements  of  Russian  workers  [Kozlov,  et  al.,  1971; 
Kravtsov,  etal.,  1971,  and  Groskaya,  et  al. ,  1972]  indicate  that  sea 
waves  approaching  a  shoreline  produce  appreciable  magnetic  fields. 
Indeed,  magnetic  fieid  magnitudes  per  unit  wave  height  appear  to  be 
anomalously  large. 

At  present,  then,  the  understanding  of  electric  and  magnetic 
effects  of  steady  ocean  currents  and  progressive  waves  provide  a  means 
of  investigating  simple  flow  fields  in  an  ocean.  Flow  fields  in  an  ocean, 
however,  are  commonly  in  a  confused  state,  and  statistical  methods 
are  required  to  describe  them  [ Longuet-Higgins,  1957].  The  energy 
spectrum  of  a  random  wave  field  then  provides  a  means  of  describ'  >g 
its  statistical  properties,  and  conversely,  measurement  of  statistical 
properties,  in  principle,  provides  a  means  of  determining  energy 
spectra. 
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Our  overall  purpose  is  twofold:  (1)  to  establish  a  correspon¬ 
dence  between  noise  spectra  of  electromagnetic  fields  induced  by  ocean 
waves  and  energy  spectra  of  both  surface  and  internal  waves,  and 
(2)  to  suggest  suitable  procedures  for  determining  wave  spectra  from 
measurements  of  electromagnetic  noise  spectra.  Correspondence 
between  electromagnetic  noise  spectra  and  ocean  wave  spectra  is  pro¬ 
vided  by  a  transfer  function,  which  is  the  electromagnetic  field  induced 
by  a  progressive  ocean  wave  of  unit  wave  height. 

Ideally,  noise  spectra  of  induced  magnetic  and/or  electric 
fields  are  simply  the  product  cf  the  squared  modulus  oi  the  transfer 
function  and  ocean  wave  energy  spectra,  because  induced  fields  are 
responses  of  a  linear  network  —  the  Earth’s  magnetic  field.  Measure¬ 
ments  must  contend  with  electromagnetic  noise  produced  by  ionospheric 
and/or  geologic  sources  as  well  as  separate  noise  spectra  produced 
by  surface  and  internal  waves.  Consequently,  properly  designed  experi¬ 
ments  lead  to  more  complicated  relations  between  spectra  and  thus 
require  rather  complete  information  on  the  transfer  function,  for  both 
design  and  interpretation. 

Our  immediate  purpose  herein,  then,  is  to  develop  a  complete 
description  of  electromagnetic  fields  induced  by  both  surface  and  inter¬ 
nal  ocean  waves.  A  subsequent  paper  presents  correspondences  between 
electromagnetic  noise  spectra  and  ocean  wave  spectra. 

We  begin  by  presenting  in  Section  II  a  general  solution  for  the 
electromagnetic  field  induced  by  an  irrotational  flow  field  using  Ampere's 
iaw  and  Faraday’s  law  and  show  that  fields  are  determined  by  an  electro¬ 
magnetic  potential  function  and  a  velocity  potential  function.  We  then 
determine  solutions  corresponding  to  progressive  waves  in  a  horizontally 
stratified  ocean. 


3 


To  develop  insight,  we  first  investigate  fields  produced  by 
progressive  waves  in  a  deep  ocean  in  Section  in  and  then  consider 
fields  produced  by  progressive  waves  in  shallow  seas  in  Section  IV. 

In  each  case,  we  present  fields  produced  by  both  surface  and  internal 
waves  in  both  polar  and  equatorial  regions.  Finally,  we  give  complete 
expressions  for  electromagnetic  potential  functions  corresponding  to 
surface  and  internal  waves ( in  the  Appendix. 


II.  ELECTROMAGNETIC  FIELDS  INDUCED 
BY  OCEAN  WAVES 


Sea  water  moving  across  the  Earth's  magnetic  field  induces 
an  electric  current  in  the  sea  that  generates  an  electromagnetic  field 
characteristic  of  the  velocity  field  in  the  sea.  Motion  of  sea  water 
with  velocity  V across  the  Earth's  constant  magnetic  field,  B,  induces 
an  electric  field,  E.  =  v"  x  B.  The  electric  current,  T,  produced 
in  the  sea  is  expressed  in  MKS  units  as 

J  =  cr(E  +  E.)  ,  (1) 

where  E  is  the  electric  field  in  a  stationary  frame  of  reference  and 
<7,  the  electrical  conductivity,  which  is  about  4  mhos/m  for  sea  water. 
Because  phase  speeds  of  ocean  waves  are  much  less  than  the  sp^cu  o  1 
light,  displacement  currents  in  the  sea  are  negligible,  and  the  mag¬ 
netic  field,  b,  of  the  induced  current  is  determined  by  Ampere's  law 
in  the  form 


5  *  b  =  V  •  (2) 

_7 

whore  uQ  =  4 nx  10  henry/m.  A  second  order  contribution  to  the 
induced  current  proportional  to  "v  x  b  is  negligible,  since  amplitudes 
of  induced  magnetic  fields  are  several  ord  rs  '  f  magnitude  smaller 
than  the  Earth’s  field.  Finally,  electric  and  magnetic  fields  arc  related 
bv  Faradv’s  law,  so  that 


(3) 


Equa'.ions  (1),  (2),  and  (3),  together  with  appropriate  boundary 
conditions,  determine  the  electromagnetic  deld,  E  and  b,  generated 
by  sea  water  motion  having  a  velocity  field  v\  Because  the  conductiv¬ 
ity  of  sea  water  is  small,  the  Lorentz  body  force,  J  x  B,  acting  on 
the  fluid  as  a  result  of  the  induced  current  is  negligibly  small  com¬ 
pared  to  pressure  and  buoyancy  forces.  Consequently,  hydromagnetic 
effects  are  absent. 

For  our  purposes,  sea  water  flow  fields  are  adequately 

described  bv  'rrotational  motion  of  an  incompressible  and  inviscid  fluid, 
so  that  the  velocity  field  is  the  gradient  of  ;  p.;  •'■entia  1  function,  v  = 

-"v  and  both  its  divergence  and  vorticity  vanish.*  The  induced  elec¬ 
tric  field,  then,  is  solenoidal  and  is  expressed  in  terms  of  the  potential 
function  as 

£.  =  -?  x  ^B  .  (4) 

In  a  sea  having  a  uniform  electrical  conduc  tivity,  the  divergence  of  E 
also  vanishes,  so  that  all  vector  fields  are  solenoidal. 

By  taking  the  curl  of  both  Equations  (2)  and  (3),  we  obtain  the 
expressions 

v  E  -  M0CT  Tt  ~~  V  *T  (5a) 

and 

-  Hoo^  =  -Mocrv  x  E.  ,  (5b) 

relatinu  the  solenoidal  field  vectors  E  and  b  to  the  induced  electric  field  E.. 

*  Vorticity  is  important  in  tidal  motions,  because  Coriolis  forces  then 
int  ku  nee  flow  fir ’ds . 

fi 


B.  HOMOGENEOUS  SOLUTION 

Because  a  solenoidal  vector  field  is  determined  by  two  scalar 
fine  tic  as  [Morse  and  Feshbach,  1953],  we  express  the  electric  field,  "E^, 
(.hot  satisfies  the  homogeneous  form  of  Equation  (5a)  in  terms  of  two 
scalar  functions,  4>j  and  and  a  unit  vector  £  directed  vertically  down¬ 
ward  in  the  form 

E^  -  v~  x  ( ^ jZ  +  z  x  v-^)  ’  (9a) 

where  both  <fj  and  satisfy  the  expression 


7 
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(9b) 


The  function  <£j  generates  a  transverse  electric  type  field,  so  called 
because  the  corresponding  electric  field  vector  is  transverse  to  the 
gradient  direction,  and  $2  generates  a  so-called  transverse  magnetic 
type  field. 

Above  the  ocean  surface,  conductivity  vanishes.  As  a  result, 
the  transverse  magnetic  component  reduces  to  an  electrostatic  type 
field;  fields  above  the  surface  are  given  bv 

*a  x*.l*  +’*.2  (lla) 

and 

db 

-q^  =  -T(z-v‘«al)  ,  (lib) 

where  both  $>aj  and  $a2  satisfy  a  Laplace  equation 

=  0  .  (11c) 

ai 
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Since  the  particular  solution  vanishes  above  the  ocean  surface,  the 
fields  E  andF  are  complete  fields, 

<1  ct 

C.  GENERAL  SOLUTION 

Vertical  electric  currents  in  the  sea  are  rapidly  extinguished 

by  accumulation  of  surface  charge.  Electric  current  is  proportional 

to  the  homogeneous  electric  field,  E^,  so  that  the  vertical  component 

of  ITh  vanishes.  Consequently,  the  potential  function  *2  vanishes,  and 

homogeneous  fields  are  transverse  electric  type  fields.  Complete 

fields  in  the  sea,  If  and  b  ,  which  are  a  sum  of  particular  and  homo- 

s  s 

geneous  fields,  are  then  given  by 

E  =?x  (♦  z  +  *B)  (12a) 

5  S 

and 

<9tT  __  __  a<t>  , 

_  =  -v(B  .  7*  +  z  •  v  *g)  +  ^Qcr-~  z  ,  (12b) 


where  the  potential  function  <t>  satisfies  Equation  (9b). 

s 

Above  the  ocean  surface,  complete  fields  are  given  by  Equa¬ 
tions  (11a)  and  (lib).  Below  the  ocean  floor,  electric  and  magnetic 
fields,  !Tc  and  b"  ,  are  expressed  as 


E  =  v  x  <t>  z 
c  c 


and 


dh 
_ c 

at 


-  "v(z 


V*  ) 
c 


+  urr 
HO  C 


dtc 

9t 


A 

Z 


(13a) 


(13b) 
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where 


0 


and  ar  is  the  conductivity  of  suboceanic  strata, 
c 


(13c) 


D.  PROGRESSIVE  WAVE  SOLUTION 

As  a  first  approximation,  we  consider  oceans  to  be  horizon¬ 
tally  stratified  and  reckon  the  vertical  coordinate,  z,  positive  with 
increasing  depth.  We  direct  the  horizontal  x-axis  northward  and  the 
y-axis  eastward,  so  that  the  Earth's  magnetic  field  vector  lies  in  the 
x-z  plane  and  is  inclined  at  an  angle  (f>  to  the  x-axis. 

Because  we  represent  oceans  as  a  composite  of  horizontal 
strata,  we  express  each  of  the  scalar  potential  functions  that  specify 
the  electromagnetic  field  as  a  Fourier  integral  of  a  potential  function 
corresponding  to  a  horizontally  progressing  ocean  wave  having  an  angu¬ 
lar  frequency  «j  and  horizontal  wave  vector  k.  We  express  the  velocity 
potential  function,  for  example,  as 

'Hz,?,  t)  =  J J*  dwdic<Kz,]T,G>)  e*  ^  -  k  r)  ^  (14) 

where  T  is  a  horizontal  coordinate  vector  The  Fourier  coefficient 
i/'(z,lc,w)  is  the  vertical  profile  of  a  wave  progressing  horizontally  in 
a  direction  k  with  an  angular  frequency  w. 

Linear  operators  relate  the  electromagnetic  field  to  potential 
functions.  In  order  to  determine  the  electromagnetic  field  produced 
by  a  specified  flow  field  in  the  sea,  then,  we  first  relate  Fourier  coef¬ 
ficients  of  the  electromagnetic  field  to  those  of  potential  functions. 

We  determine  the  relation  by  finding  the  electromagnetic  field  produced 
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by  a  progressive  wave  in  which  potential  functions  have  the  form  of 
the  integrand  of  Equation  (14). 

For  a  progressive  wave,  the  profile  of  the  velocity  field, 
v(z,lT,«),  is  a  reduced  gradient  of  the  profile,  »Hz,F,ti>),  of  its  poten¬ 
tial  function 

v  =  ,  (15a) 


where  the  reduced  gradient  operator,  v^,  is  defined  by 


a  d 
z  dz 


-  ik 


(15b) 


Since  the  velocity  potential  function  satisfies  a  Laplace  equation,  its 
profile  in  a  progressive  wave  satisfies  the  Helmholtz  equation 

^4  -  k2^  =  0  .  (15c) 

dzz 

To  facilitate  further  calculation  and  interpretation,  we  express 
7  as  a  sum  of  circularly  polarized  components  in  the  form 


v 


~p=r-  [(1  +  p)0  -  (1  -  p)0*j 
v2 


(16a) 


where 


t?  \  _  1  d  In 

p(z,  K, <y)  -  -  £ 

v  =  -p=~  (z  +  ifi)  , 
/2 


(16b) 

(16c) 
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k  is  a  unit  vector  in  the  direction  of  wave  propagation,  k  is  a  wave 
number,  and  v*  is  the  complex  conjugate  of  0.  The  unit  vector  u  has 
the  convenient  property  that  Ox  n  =  if),  where  n  =  $  x  £,  so  that  v  n=0. 
Circular  polarization  of  t)  is  left-handed  with  respect  to  n,  and  that  of 
its  complex  conjugate  is  right-handed. 

The  velocity  vector  of  a  progressive  wave,  then,  is  elliptically 
polarized  in  a  vertical  plane  containing  the  wave  vector.  Eccentricity 
of  the  ellipse  and  sense  of  rotation  with  respect  to  the  direction  of  the 
unit  vector  n,  normal  to  the  plane  of  polarization,  are  determined  by 
the  value  of  the  function  p(z,"E,co),  which  we  call  the  polarization  coef¬ 
ficient. 

From  Expression  (12a)  that  relates  the  electric  field  in  the 
sea  to  a  velocity  potential  and  a  potential  function  *  ,  satisfying  Equa- 
tion  (9b),  we  find  that  the  electric  field  profile,  T  (z,lcja>),  of  a  pro- 

9 

gressive  wave  is  given  by  the  expression 


(17a) 


where  <f>  (z,k,f>)  is  the  progressive  wave  profile  of  the  potential  func- 
s 

tion  <t>  and  so  satisfies  the  Helmholtz  equation 
s 


(17b) 


2  2  -► 
with  a  =  k  +  i^cro.  By  expressing  the  vector  product  yx  £T  as  a 

sum  of  components  along  and  transverse  to  the  unit  vector  n,  we 

write  Equation  (17a)  as  a  sum  of  transverse  electric  and  electrostatic 

type  fields  in  the  form 
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(18a) 


7  =  [Ik*  +  0- (? x  n)]n  -  (B-nXiJxS)  . 

s  s 

The  profile  of  the  electric  current  produced  in  the  sea  by  a  progressive 
wave  is  simply 

Ta  =  ik<r<^n  .  (18b) 

By  taking  the  curl  of  the  electric  field,  we  obtain  the  expression 

i (oJs  =  n  xvr  [ik<£s  +!*(17xn)]  ,  (19) 

for  the  profile  of  the  magnetic  field. 

By  writing  the  reduced  gradient  of  6  in  terms  of  circularly 

S 

polarized  components,  we  express  magnetic  and  electric  field  profiles 
in  terms  of  circularly  polarized  components  as 


=  ik 


(f> 

—  (1  +  p  )  +  *(  1  +  p)(B-0)| 

VT 


and 


+  ik 


(1  -  P  )  +  <H1  -  p)(B*  P*) 

VT  8 


(20a) 


ikU  +  —  [(1  +  p)(B*i))  +  (1  -  p)(B*  t)*)]!n 

S  y[2  ) 


-  —  (B.n)[(l  +p)£  +  (1  -p)0*]  .  (20b) 

VT 
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We  observe  from  Equation  (20b)  that  the  electrostatic  part  of  the  elec¬ 
tric  field  vanishes  in  polar  regions  where  the  Earth's  magnetic  field  is 
nearly  vertical  so  that  15* ri  vanishes.  The  field  is  then  a  pure  trans¬ 
verse  electric  type  field.  In  equatorial  regions,  the  field  is  a  pure 

electrostatic  type  field  ?*en  5  =  Bn,  so  that *5* D  and  hence  6  vanish. 

s 

We  also  note  that  the  ratio  of  the  transverse  component  of  the  electric 
field  to  the  vertical  component  of  the  magnetic  field  is  the  phase  speed 
of  an  ocean  wave;  namely, 


(20c) 


In  conducting  strata  below  the  oceans,  the  velocity  potential 

vanishes,  and  electric  and  magnetic  fields  are  transverse  electric 

type  fields.  Their  profiles  are  expressed  by  Equations  (20b)  and  (20a) 

with  4>=  0  and  <t>  replaced  by  a  potential  profile  $  .  The  profile  <f> 
s  c  c 

satisfies  a  Helmholtz  equation  like  Equation  (17b)  with  a  constant, 

2 

a  ,  corresponding  to  the  conductivity  of  suboceanic  strata. 

V 

Above  the  ocean  surface,  the  electric  and  magnetic  field  pro¬ 
files,  and  /f ,  are  given  in  terms  of  potential  profiles  and  <f>a2 
by  the  expressions 


e 

a 


ik*aia  +  Va2 


(21a) 


and 


(21b) 
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irpnn 


where  both  4>aj  and  <6^  satisfy  the  equation 

dV  2 

-  kV  =  0  .  (21c) 

dz 1  ai 

Because  fields  must  vanish  far  above  the  surface,  potential  profiles 
are  of  exponential  form,  so  that 


T  <f>  .  =  k  0  .  v*  .  (21d) 

r  ai  aiv 

Field  profiles  above  the  surface  are  expressed  in  terms  of  circularly 
polarized  components  as 

=  +  ^2"  k^a2  v*  (22a) 


and 


(22b) 


Whatever  the  polarization  of  the  velocity  field,  the  magnetic  field 
produced  above  the  surface  by  a  progressive  wave  is  circularly 
polarized  with  a  right-handed  sense  of  rotation. 

In  each  horizontal  layer  (air,  strata  within  a  sea,  and  sub- 
oceanic  strata),  specification  of  two  arbitrary  constants  determines 
the  electromagnetic  field  produced  by  a  progressive  wave.  Requiring 
continuity  of  horizontal  components  of  electric  and  magnetic  field  vec¬ 
tors  at  interfaces  together  with  requiring  fields  to  vanish  far  enough 
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above  and  below  an  ocean  specify  the  two  constants  in  each  layer. 
Further  description  of  the  field,  then,  requires  specifying  a  model  to 
represent  an  ocean.  In  order  to  develop  insight,  we  begin  with  the 
simplest  stratification  mode!  —  a  homogeneous,  deep  ocean  —  and  pro¬ 
ceed  to  multilayered  models  representing  shallow,  stratified  seas. 


in.  ELECTROMAGNETIC  FIELDS  INDUCED  BY 
PROGRESSIVE  WAVES  IN  A  DEEP  OCEAN 


In  this  section,  we  describe  electromagnetic  fields  produced 
by  progressive  waves  in  an  ocean  that  is  several  wave  lengths  deep 
(kD»  1).  We  first  consider  surface  waves  on  an  homogeneous  ocean 
and  then  internal  waves  propagating  along  a  sharp  thermocline.  In 
each  case,  we  suppose  that  wave  amplitudes  are  much  smaller  than 
a  wavelength  and  so  restrict  consideration  to  a  linear  regime. 


A.  SURFACE  WAVES 


In  a  deep  ocean,  the  velocity  potential  profile,  and  the  electro- 
magnetic  potential  profile  <i>3  decrease  exponentially  with  depth  as  e 

m(XZ 

and  e  ,  respectively,  since  all  fields  vanish  at  great  depths.  The 

polarization  coefficient  of  the  velocity  profile  is  unity,  and  p  =  a/k, 

s 

where 


(23a) 


with 


2 

5 


UqCTW 


(23b) 


The  quanity  b/y/2  is  the  ratio  of  wavelength  of  an  ocean  wave  to  wave- 

1  /2 

length  of  an  electromagnetic  wave  in  the  sea,  \  ,  where  \  =  2r(2/u  era) 

v  v  VJ 

[Kraichman,  1970]. 


From  Equations  (20a)  and  (20b),  we  then  find  that  electric  and 
magnetic  field  profiles  in  a  deep  ocean  are  given  by  the  expressions 


?  =  ik  [*oe-“z  +  *0>/r( B-0)  e-kz]B 


-  Utifi  (B.ft)  e"kz  yJT  0 


(24a) 


and 


m 


& 


H 

I  *•(«-■) 


«  A\  -kz 


e  +  2+JT  (B-O)e 


-az 

e  0* 


(24b) 


where  $  and  are  values  of  <l»  and  at  the  surface,  z  =  0.  The  pro- 
0  0  s  r 


file  of  the  current  induced  in  the  sea,  j  ,  is  expressed  as 

s 


•m  •  i  i  -az  a 

js  =  <nk0Q  e  n 


(24c) 


Above  the  surface  (z  <  0),  field  profiles  are  expressed  in  terms 

of  two  constants,  and  A, ,  as 

o  1 


=  ik  e^z  (AQn  -  iA^y[Tv*) 


(25a) 


and 


(iH 


ikA  e^z  yj%"  0* 


(25b) 
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By  imposing  continuity  of  horizontal  field  components  at  the 
surface,  we  obtain  the  expressions 


■  -(«Te)  *„*<*•& 

Ao  ■  (jrlk^S-O) 


and 

Aj  =  i  ^  (B-rt) 


(26a) 


(26b) 


(26c) 


relating  the  constants  <t>Q,  Aq,  and  Aj  to  ^Q.  For  ocean  wave  ampli¬ 
tudes  that  are  much  less  than  wavelengths,  the  profile  of  wave  ampli¬ 
tude,  £,  is  expressed  in  terms  of  the  velocity  potential  profile  by  the 
kinematical  relation 

io>£  =  kp^  ,  (27a) 

so  that  wave  height  at  the  surface,  f  ,  of  a  progressive  wave  in  a  deep 
ocean  is  given  in  terms  of  by  the  expression 

i<u£o  =  .  (27b) 

Fields  in  the  sea,  then  are  proportional  to  surface  wave 
heights  and  are  expressed  by 
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i  t  uni  ■!»■■■  wjj u wmtmm SPWWff 


?s  ■ 


■  (^Tk) 


•(a-k)z 


(B*$)n 


and 


-  (B-n)u! 


j3s=  -kfoe_kz2(g.C)| 


1  -  e 


-(a  -  k)z 


-  k)z  Aj 


(28a) 


(28b) 


below  the  surface  (z  >  0)  and  by 


ra  =-"f0  +  (B-fi)  0* 


and 


A  = -k£ekz  2(B*0) 

d  U 


y 


(28c) 


(28d) 


above  the  surface  (z  <  0).  The  profile  of  the  current  induced  in  the  sea 
is  given  by 


-r  „  -kz 

i  =  (hoc  e 
Js  o 


(s?r)  e'(“ '  kV  •  <»> 
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Bernoulli's  equation  in  linearized  form  [Lamb,  1932]  gives 


the  relation 


2 

iwf  =  Q-  <|/  . 

o  g  ro 


(30a) 


between  wave  height  and  the  velocity  potential  profile  at  a  free  surface, 

2 

where  g  =  9.8  m/sec  .  The  kinematical  relation  given  by  Equation  (27b) 

2 

together  with  Equation  (30a),  then,  give  the  dispersion  law,  co  =  gk, 
for  surface  waves  in  a  deep  ocean.  By  using  the  dispersion  law,  we 
express  the  square  of  the  ratio  of  ocean  wavelength  to  electromagnetic 
wavelength  as 


2k 


so  that 


(5xl0"7) 


3/2 


(30b) 


(30c) 


where  \  is  wavelength  measured  in  meters.  For  ocean  wavelengths 
less  than  a  few  kilometers,  corresponding  to  wave  periods  less  than 
a  few  tens  of  seconds,  wavelengths  of  ocean  waves  in  a  deep  ocean  are 
much  less  than  thos°  of  electromagnetic  waves.*  The  ratio  a/k,  deter¬ 
mined  by  Equation  (23a),  is  then  given  by 


*  In  an  ocean  of  depth  D,  phase  speeds  of  surface  waves  aie  less 
than  i/gDfLamb,  1932],  so  that  8^/2  <  4 X  v/gP  x  10"^.  Because 
ocean  depths  are  of  the  order  of  10  km  at  most,  sV 2  <  X  x  10“^. 


21 


a  ~  i  ,iP 

r  «  ■*•  +  l  “7T 


(30d) 


to  a  good  approximation. 


Using  the  approximate  form  for  the  ratio  a/k,  we  obtain  first 


order  expressions  in  8  for  fields  in  and  above  the  ocean;  namely, 


T  =  e"kz 

S  0 


v»[' 


i-^-(l  +  2kz)(B*p)fi  -  (B'n)i) 


Ola) 


-  -ikf  e'kz  4-(B-0)  (2kz0  +  0*)  , 


(31b) 


=  -  o,(o  ekzy/t  [i  (B-0)n  +  (B-n)  £*]  ,  (31c) 


and 


-  ik£Q  e^z  -y  (B-v)  v* 


(31d) 


The  induced  current  profile  in  the  sea  is  expressed  to  first  order  as 


ctoj £  e"kz  [1  -  i  4r  (1  +  2kz)]^2  (B*i))n  .  (31e) 


We  note  that  the  transverse  electric  component  of  the  electric 

2 

field  and  the  magnetic  field  vanish  as  S  becomes  negligible,  so  that  to 
zero  order  the  field  is  an  electrostatic-like  field  produced  by  surface 
charge.  The  zero  order  current  is  driven  by  the  transverse  component 
of  the  electric  field  induced  by  sea  water  motion.  In  polar  regions,  how¬ 
ever,  B-nsO,  so  that  surface  charge  and  hence  the  zero  order  electric 
field  vanish,  and  only  the  transverse-electric  type  field  remains. 
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From  Equation  (31b),  we  find  that  magnitudes  of  induced  mag¬ 
netic  fields,  to  first  order,  are  proportional  to  phase  speeds  of  surface 
waves,  since 


(32a) 


Similarly,  magnitudes  of  the  transverse-electric  components  of  elec¬ 
tric  fields  are  proportional  to  the  square  of  phase  speeds.  Because 
phase  speeds  of  surface  waves  on  a  deep  ocean  are  inversely  propor¬ 
tional  to  frequency  or,  equivalently,  to  the  square  root  of  wave  num¬ 
ber,  field  magnitudes  increase  with  decreasing  frequency  or  wave 
number . 

We  remark  that  an  increasing  importance  of  magnetic  induc¬ 
tion  effects  with  decreasing  frequency  is  contrary  to  customary  notions. 
Ordinarily,  induction  effects  are  more  pronounced  at  higher  frequencies. 
An  increase  in  magnitude  of  induced  fields  with  decreasing  frequency 
does  not  continue  indefinitely,  however,  because  at  zero  frequency  the 
sea  in  quiescent  and  induced  fields  vanish. 

To  illustrate  the  dependence  of  field  magnitudes  on  wave  num¬ 
ber,  we  consider  fields  induced  at  the  surface  in  polar  regions.  For 
ocean  wavelengths  that  are  much  smaller  than  electromagnetic  wave¬ 
lengths,  we  find  from  Equation  (31  q)  that  the  magnetic  field  strength 
induced  at  the  surface  by  a  one  meter  high  wave  in  polar  regions  is 
given  by 


K\ 


(32b) 


and  so  increases  with  decreasing  wave  number.  For  ocean  wavelengths 


that  are  much  greater  than  electromagnetic  wavelengths,  however,  we 
find  from  Equation  (28d)  that  magnetic  field  strength  at  the  surface  is 
proportional  to  wave  number;  namely 


Bk  y/I 


(32c) 


and  so  decreases  with  decreasing  wave  number.  In  each  case,  cor¬ 
responding  electric  field  magnitudes  are  determined  by  the  relation 


|e0|  =  |bj  viTat  . 


(32d) 


Dependence  of  magnetic  and  electric  field  magnitudes  on  wave 
number  is  delineated  in  Figure  1  by  the  heavy  solid  and  dashed  curves, 
respectively.  The  light  dashed  lines  are  asymptotes  determined  by 
Equations  (32b)  and  (32c).  Curves  in  Figure  1  correspond  to  field  mag¬ 
nitudes  induced  at  the  surface  by  a  one  meter  high  wave  in  polar  regions, 
where  the  Earth’s  field  strength  is  about  6. 24  x  10  gammas  [Chapman, 
1962].  As  is  evident  in  Figure  1,  field  magnitudes  attain  maximum 
values  at  wavelengths  near  10  km,  where  electromagnetic  and  ocean 
wavelengths  are  comparable.  The  magnitude  of  the  magnetic  field 

attains  a  maximum  value  of  about  10  gamma/m,  and  the  magnitude  of 

2 

the  electric  field,  a  value  of  about  1  microvoit/m  . 

Nonetheless,  ocean  wavelengths  of  several  kilometers  or 
greater  are  comparable  to  or  exceed  ocean  depths,  so  that  a  deep 
ocean  approximation  no  longer  applies,  and  phase  speeds  of  surface 
waves  approach  the  limiting  value  y/gD.  Determination  of  fiela  mag¬ 
nitudes  induced  by  surface  waves  with  wavelengths  exceeding  a  few 
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kilometers,  corresponding  to  wave  periods  exceeding  a  few  tens  of 
seconds,  requires  accounting  for  wave  reflection  from  suboceanic 
strata  —  a  topic  that  we  consider  in  Section  IV. 

When  a  deep  ocean  approximation  is  adequate,  however,  ocean 
wavelengths  are  much  less  than  electromagnetic  wavelengths,  so  that 
Equations  (31a)  to  (31d)  adequately  describe  field  profiles  in  a  deep 
ocean.  For  purposes  of  discussion,  we  first  consider  polar  regions. 
The  magnetic  field  vector  above  the  surface  in  polar  regions  is  cir¬ 
cularly  polarized  iu  a  right-handed  sense  with  respect  to  n,  and  its 
magnitude  decreases  exponentially  with  altitude  in  accordance  with  the 
relation 


ibl  =  |b0|  ekZ  ,  (33a) 

with  |bQ|  determined  by  Equation  (32b). 

Below  the  surface,  the  magnetic  field  vector  is  ellipitically 
polarized  with  the  major  axis  of  the  ellipse  aligned  vertically.  The 
sense  of  rotation  is  right-handed  at  depths  less  than  x/4tt.  At  a  depth 
of  a/47 r,  the  field  is  linearly  polarized  in  a  vertical  direction.  Below 
a  depth  of  \/47r,  the  field  is  again  elliptically  polarized  with  the  major 
axis  of  the  ellipse  aligned  vertically,  but  the  sense  of  rotation  is  left- 
handed.  At  great  depths  polarization  becomes  circular  with  a  left- 
handed  sense  of  rotation.  The  magnitude  of  the  magnetic  field  below 
the  surface  decreases  with  depth  according  to  the  relation 

Ibj  =  jbj  e  ,  (33b) 
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Profiles  of  magnetic  field  strength  and  polarization  induced 
in  polar  regions  by  a  one  meter  high  wave  having  a  wavelength  of  10  m 

are  shown  in  Figure  2.  Altitude  and  depth  are  measured  in  units  of  a 
wavelength. 

In  polar  regions,  surface  charge  vanishes,  and  electric  fields 
are  transverse  electric  type  fields.  The  magnitude  of  the  electric  field 
decreases  above  the  surface  as 

®  ■  l®ol  ekz  (34a) 

and  below  the  surface  as 


®  *  IEJ  (1  +  2k  z)  e*kz  , 

where  the  magnitude  at  the  surface,  JEJ,  is  given  by 


(34b) 


(34c) 


The  magnitude  of  the  electric  current  density  in  polar  region 
decreases  exponentially  with  depth  in  accordance  with  the  relation 


iJj  -  cr(o^QB  e 


(34d) 


Magnitudes  of  current  densities  at  the  surface  are  inversely  propor¬ 
tional  to  wave  period  and  for  a  one  meter  high  wave  in  polar  regions, 
range  from  about  1,5  milliamps/m2  at  a  wave  period  of  1  sec  to  0.015 
milliamperes/m2  at  a  period  of  100  seconds. 
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In  mid-latitude  and  equatorial  regions,  a  circularly  polarized 
electrostatic  component,  resulting  from  surface  charge,  is  added  to 
the  electric  field.  The  sense  of  rotation  of  the  electrostatic  component 
is  right-handed  with  respect  to  n  above  the  surface  and  left-handed 
below  the  surface.  Its  magnitude  decreases  exponentially  both  above 
and  below  the  surface  as 

If  I  =  e’k^  (35a) 

V  w 

where 

iB-ftl  =  Bpcoa^inl  V»+3sin2«.  ,  (35b) 

B  is  the  Earth’s  field  strength  in  polar  regions;  0,  the  angle  between 
1? 

the  direction  of  wave  propagation  and  magnetic  north;  <t>,  magnetic 
latitude,  and  tan  $  =  2  tan$>  relates  dip  angle,  <f>,  and  latitude.  A  one 
meter  high  wave  propagating  across  magnetic  meridians  at  the  mag¬ 
netic  equator  induces  an  electrostatic  field  strength  at  the  surface  of 
0. 28  millivolts/m  for  a  wave  period  of  1  sec  and  0.028  millivolts/m 
for  a  wave  period  of  10  sec. 

In  equatorial  regions,  magnitudes  of  the  transverse  electric 
component  of  the  electric  field  and  of  the  magnetic  field  are  less  than 
their  magnitudes  in  polar  regions  by  a  factor  of  cos  0/2.  At  inter¬ 
mediate  latitudes,  magnitudes  of  surface  fields  are  less  than  magni¬ 
tudes  in  polar  regions  by  a  factor  of 

1/2 

R  =  — °-y-  [(1  +  3  sin^  <$(tan^  <f>  +  cos^  d)]  .  (35c) 

Polarization  of  the  magnetic  field  vector,  however,  is  the  same  as  in 
polar  regions. 
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B.  INTERNAL  WA¥ES 

For  a  deep  ocean  containing  a  sharp  thermocline  at  a  depth  d, 
we  express  the  profile  of  the  velocity  potential  for  a  progressive  wave 
as 

^1  =  *o  (cosk  kz  “  ^  sinh  kz)  ,  (36a) 

above  the  thermocline  (0<z<d),  and  as 

*  =  ^  e“k(z“d)  (36b) 

i  d 

below  the  thermocline  (z>d),  where  *PQ  and  ^  are  values  of  respective 
profiles  at  the  surface  (z  =  0)  and  at  the  thermocline  interface  (z  =  d). 
Bernoulli’s  equation  at  a  free  surface,  as  given  by  Equation  (30a), 
together  with  the  kinematical  relation  expressed  by  Equation  (27a)  give 
the  relation 

i)2  =  &j2/gk  (37) 


for  the  value  of  the  polarization  coefficient  at  the  surface. 

2 

For  surface  waves,  17  is  unity,  corresponding  to  the  disper- 

2  2 
sion  law  <u  =  gk.  For  internal  waves,  however,  17  is  determined  by 

the  relation 


Wp2 

Pj/Pg  +  coth'kd 


(38a) 


which  assures  continuity  of  pressure  and  the  vertical  component  of 
velocity  at  a  thermocline  interface.  Sea  water  density  below  a  thermo¬ 
cline,  Pg,  exceeds  that  above  a  thermocline,  Pj,  by  a  small  amount 


Ap.  To  a  good  approximation, 


-3 

where  A  p/p  is  of  the  order  of  10  [Phillips,  1969].  Requiring  contin¬ 
uity  of  the  vertical  component  of  velocity  alone  at  a  thermocline  gives 
the  relation 


$  sinh  kd  (1  -  tf  coth  kd) 


(38c) 


between  values  of  respective  profiles  at  the  surface  and  thermocline 
interface .  The  horizontal  component  of  velocity  is  discontinuous  at  a 
sharp  thermocline. 

We  express  the  profile  4>s  of  the  electromagnetic  potential  in 
terms  of  two  constants,  and  Cg,  as  the  sum  of  exponentials 

+gl  *  Cj  e“z  +  C2  e'az  ,  (39a) 

above  the  thermocline  (0<z<d),  and  in  terms  of  a  constant  Cg  as 

*2  =  C3e'“(z'd)  (39b) 

below  the  thermocline  (z>d),  since  </>  vanishes  at  great  depths.  From 

s 

Equations  (20a)  and  (20b),  we  then  obtain  the  expressions 


29 


and 


'  jrh0’*"-  +  c2(1+^e'"z 


+  '/'o(l+7?2V5”(B*0)e"kzp 
+  *0(l~A^(B>)ekz  U*  ,  (40b) 


for  electric  and  magnetic  field  profiles  above  the  thermocline,  and 
the  expressions 

%2  =  *  [*„  +  V5  (B-0)e-k<2-d)J  ft 

-  ik^T(B-n)e"k(z"d)  v  (40c) 

and 
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(tfs2  ♦  2Vr<ff-0,e-k<~»}> 

+ j=  [s'1- ^iz'd)y  .  («« 

for  electric  and  magnetic  field  profiler?  below  the  thermocline.  As 
before,  field  profiles  above  the  surface  are  expressed  in  terms  of  two 
constants,  AQ  and  Aj,  by  Equations  (25a)  and  (25b). 

By  imposing  continuity  of  horizontal  components  of  fields 

both  at  the  surface  and  at  the  thermocline,  we  express  the  constants 

o 

Cp  C2,  Cg,  Aq,  and  Aj  in  terms  of  ^  and  "0  .  Because  the  resulting 

expressions  are  cumbersome,  however,  we  list  them  in  the  Appendix. 

Here,  we  consider  two  approximations  that  allow  considerable  simpli- 

o 

cation  in  expressing  field  profiles.  First,  for  internal  waves,  V  «1, 

9 

so  we  everywhere  neglect  as  compared  to  unity.  As  a  consequence, 
wave  height  at  the  free  surface  is  negligibly  small,  and  the  potential 
function  values  ^  and  ^  are  related  to  wave  height  at  the  thermocline, 
or  thermocline  displacement  by  the  expression 

=  S  "  *0  sinh  kd  ’  (41) 

which  we  obtain  from  the  kinematical  relation  and  Equation  (38c) . 

Second,  the  ratio  of  wavelength  of  an  internal  wave  to  wave¬ 
length  of  an  electromagnetic  ’vave  is  small  in  a  deep  ocean,  as  for  sur¬ 
face  waves.  We  express  the  square  of  the  ratio  in  terms  of  phase  speed 
of  an  internal  wave  as 
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(42a) 


/»  -  ¥  (I)  ■ 

and,  from  Equation  (38b),  we  express  the  phase  speed  as 


Thus,  internal  wave  phase  speeds  are  less  than  v47a  p/p)d  ,  and  so 

82/ 2  <  (4  x  10‘7)\yg(AP/P)d  ,  (42c) 


where  X  is  wavelength  of  an  internal  wave.  Since  thermocline  depth 

-3 

cannot  exceed  ocean  depths  (d<10  km)  and  A  P/p  is  of  the  order  of  10  , 
82/2  <  (4xlO‘6)A  ,  (42d) 


so  that  the  ratio  is  small  for  wavelengths  less  than  a  few  tens  of  kilo¬ 
meters.  For  a  deep  ocean  approximation,  wavelengths  are  always  of 
the  order  of  or  less  than  a  few  kilometers.  Wavelengths  of  internal 
waves,  then,  are  much  smaller  than  wavelengths  of  electromagnetic 
waves  in  a  deep  ocean,  and  so  to  a  good  approximation 

«/k  »  l+i82/2  ,  (42e) 


as  for  surface  waves. 

2 

By  neglecting  v  compared  to  unity  and  using  the  first  order 

approximation  for  <*/k,  we  express  electric  and  magnetic  field  profiles 
2 

to  first  order  in  8  as 
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*a  "  (f)<K>  S 


(43a) 


and 


\  ■  “4  t(sj srar) ekz 


(43b) 


above  the  surface  (zSO);  as 


rsi  ■  t)<Ki>» 


2  sinh  kd 


t-xtt  (B-ft)v^_(e'kzS>  +  ekz£»  (43c) 


and 


2  vl 

jg  =  ik(  5  / sinh  kz \ 

psl  m?d  2  \sinh  kd  / 


(B*0*)  +  (5.0)  5 


Ax  kz  e 


-kz 


sinh  kz 


A 

V 


kz 


A* 

V* 


(43d) 


above  the  thermocline  (0iz£d);  and  as 


rs2  =  (k^'  ?2)fi  +  '»fde-k(z;‘d) 


(43e) 


and 


r  a  ikf  i_  e“k(z«d) 

^s2  2  6 


/.  ,  -kd 

®4°*>  +(mr 


-  2k  (z-d|(B* 


(B-0)0* 


(43f) 


below  the  thermocline  (z>d). 
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For  purposes  of  discussion,  we  first  consider  polar  regions. 
From  Equation  (43b),  we  find  that  the  magnitude  of  the  magnetic  field 
produced  at  the  surface  by  a  thermocline  displacement  of  one  meter  in 
polar  regions  is  given  by 


The  corresponding  magnitude  of  the  electric  field  at  the  surface  is 
given  by 


which  we  obtain  from  Equations  (43f)  and  (43e) . 


(44d) 
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Magnitudes  of  the  magnetic  field  induced  at  the  surface  and 
at  a  100  meter  deep  thermocline  by  a  thermocline  displacement  of  one 
meter  in  polar  regions  are  shown  as  functions  of  kd  by  the  solid  curves 
in  Figure  3.  Corresponding  electric  field  magnitudes  are  delineated 
by  dashed  curves.  We  observe  that  field  magnitudes  increase  as  wave 
number  decreases,  as  for  surface  waves.  An  eventual  decrease  in 
field  magnitudes  at  still  smaller  wave  numbers  is  not  shown  in  Fig¬ 
ure  3,  as  internal  wavelengths  then  become  comparable  to  or  exceed 
ocean  depths  and  a  deep  ocean  approximation  is  no  longer  appropriate.  * 

For  shallow  thermoclines  (kd  <  0.1)  in  a  deep  ocean,  field 
magnitudes  are  effectively  independent  of  wave  number,  as  is  evident 
in  Figure  3.  And  magnetic  field  magnitudes  increase  approximately 
as  the  square  root  of  thermocline  depth,  as  is  evident  from  Equation  (44a), 
because  phase  speeds  of  internal  waves  increase  in  proportion  to  the 
square  root  of  thermocline  depth  for  shallow  thermoclines.  Electric 
field  magnitudes  increase  approximately  in  proportion  to  thermocline 
depth.  Magnitude  of  the  transverse  electric  field  component  is  quite 
small,  however,  because  phase  speeds  of  internal  waves  are  of  the 
order  of  1  m/sec  or  less. 

For  deep  thermoclines  (kd  >  1),  field  magnitudes  at  the  sur¬ 
face  decrease  rapidly  with  increasing  wave  number,  so  that  fields 
induced  by  internal  waves  having  wavelengths  of  the  order  of  or  less 
than  a  thermocline  depth  are  small  at  the  surface.  Sea  water  above 
a  thermocline,  in  effect,  acts  as  a  low  pass  filter.  Nonetheless,  we 
note  that  electric  field  magnitudes  at  the  surface  exceed  those  at  the 
thermocline  interface.  Electric  field  strength  increases  above  the 
thermocline  because  of  changing  polarization  of  the  magnetic  field. 


*  We  consider  wave  reflection  from  suboceanic  strata  in  Section  IV. 


Above  the  surface,  the  magnetic  field  vector  is  always  cir¬ 
cularly  polarized  in  a  right-handed  sense  with  respect  to  n,  and  its 
magnitude  decreases  exponentially  with  altitude.  Below  the  surface, 
magnetic  field  profiles  differ  markedly  for  shallow  and  deep  thermo- 
clines. 

Profiles  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  by  a  thermocline  displacement  of  one  meter  in  polar  regions 
are  shown  in  Figure  4  for  a  shallow  thermocline  corresponding  to 
kd  =  0.1  and  a  thermocline  depth  of  100  m.  Depth  is  measured  in 
units  of  a  thermocline  depth. 

We  observe  from  the  polarization  diagram  in  Figure  4  that 
one  polarization  reversal  occurs  above  the  thermocline  and  two  below. 
Polarization  begins  as  circular  with  a  right-handed  sense  of  rotation 
at  the  surface,  becomes  linear  with  a  vertical  alignment  at  a  depth  of 
d/4,  again  becomes  circular  but  with  a  left-handed  sense  of  rotation 
at  a  depth  of  d/2,  and  is  elliptical  with  a  horizontal  elongation  and  a 
left-handed  sense  of  rotation  at  the  thermocline  depth,  d.  Below  the 
thermocline,  horizontal  elongation  of  the  polarization  ellipse  continues 
until  the  field  becomes  linearly  polarized  with  a  horizontal  alignment 
at  a  depth  of  5.5  d.  Thereafter,  the  sense  of  rotation  again  becomes 
right-handed,  and  horizontal  elongation  of  the  polarization  ellipse 
lessens  until  polarization  once  again  becomes  circular  at  a  depth  of 
10.5  d.  Vertical  elongation  then  grows  until  polarization  becomes 
linear  with  a  vertical  alignment  at  a  depth  of  15.5  d.  Thereafter,  the 
sense  of  rotation  becomes  left-handed,  elongation  rapidly  lessens,  and 
so  polarization  is  circular  with  a  left-handed  sense  of  rotation  at  great 
depths. 
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Because  of  frequent  polarization  changes  with  increasing  depth, 
profiles  of  vertical  and  horizontal  components  of  the  magnetic  field 
exhibit  maxima  and  vanish  at  depths  corresponding  to  appropriate 
states  of  linear  polarization.  We  note,  however,  that  the  vertical  com¬ 
ponent  is  constant  above  the  thermocline,  so  that  the  electric  field  mag¬ 
nitude  is  also  constant. 

Profiles  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  by  a  thermocline  displacement  of  one  meter  in  polar  regions  are 
shown  in  Figure  5  for  a  deep  thermocline  corresponding  to  kd=4  and  a 
thermocline  depth  of  100  m.  We  observe  from  the  polarization  diagram 
in  Figure  5  that  one  polarization  reversal  occurs  at  the  thermocline, 
one  above  the  thermocline,  and  one  below  the  thermocline.  Because 
the  magnetic  field  is  linearly  polarized  with  a  horizontal  alignment  at 
the  thermocline  and  circularly  polarized  at  the  surface,  the  vertical 
component  of  the  field  vanishes  at  the  thermocline  and  increases  above 
the  thermocline.  As  a  result,  the  electric  field  strength  at  the  sur¬ 
face  exceeds  that  at  the  thermocline,  as  we  noted  previously. 

In  equatorial  regions,  magnitudes  of  the  magnetic  field  and 
transverse  component  of  the  electric  field  at  the  surface  are  less  than 
surface  magnitudes  in  polar  regions  by  a  factor  of  cos  d/2,  where  0 
is  the  angle  between  the  direction  of  internal  wave  propagation  and 
magnetic  north.  Furthermore,  the  electric  field  contains  an  electro¬ 
static  component  that  is  linearly  polarized  at  the  surface  with  a  ver¬ 
tical  alignment  and  circularly  polarized  at  and  below  the  thermocline, 
with  a  left-handed  sense  of  rotation.  Magnitudes  of  the  electrostatic 
components  are  proportional  to  sin  0. 

Magnetic  field  magnitudes  produced  at  the  surface  and  at  a 
100  meter  deep  thermocline  by  a  thermocline  displacement  of  one  meter 


propagating  along  magnetic  meridians  (0=0)  in  equatorial  regions  are 
shown  as  functions  of  kd  by  the  solid  curves  in  Figure  6.  Magnitudes 
of  the  electrostatic  component  of  the  electric  field  corresponding  to 
propagation  across  magnetic  meridians  (0=W 2)  in  equatorial  regions 
are  delineated  by  dashed  curves  in  Figure  6. 

We  observe  that  magnetic  field  strengths  at  the  surface  and 
thermocline  interface  are  equal  for  shallow  thermoclines  in  equatorial 
regions.  Magnitude  of  the  electrostatic  field  at  the  surface  decreases 
as  kd  increases,  but  magnitude  of  the  electrostatic  field  at  the  thermo¬ 
cline  interface  increases  as  kd  increases. 

Profiles  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  by  a  thermocline  displacement  of  one  meter  propagating  along 
magnetic  meridians  in  equatorial  regions  are  shown  in  Figure  7a  for 
a  shallow  thermocline,  corresponding  to  kd=0.1  and  a  thermocline 
depth  of  100  m.  Figure  7b  illustrates  magnitude  and  polarization  pro¬ 
files  for  the  electrostatic  field  corresponding  to  propagation  across 
magnetic  meridians.  For  a  shallow  thermocline  we  note  that  field 
magnitudes  are  constant  above  the  thermocline.  The  electrostatic 
field,  however,  vanishes  above  the  surface  because  internal  waves 
produce  negligible  surface  displacements.  Polarization  of  the  mag¬ 
netic  field  is  like  that  of  a  surface  wave,  as  is  that  of  the  electrostatic 
field  below  the  thermocline. 

Magnitude  and  polarization  profiles  of  magnetic  and  electro¬ 
static  fields  corresponding  to  a  deep  thermocline,  kd=4,  are  shown  in 
Figures  8a  and  8b,  respectively.  For  a  deep  thermocline,  we  note 
that  magnetic  field  profiles  are  nearly  symmetrical  with  respect  to 
the  thermocline  interface.  Because  the  electrostatic  field  vanishes 
above  the  surface,  its  polarization  above  the  thermocline  becomes 
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linear  with  a  vertical  alignment  at  surface.  The  polarization  below 
the  thermocline,  however,  is  circular  with  a  left-handed  sense  of 
rotation. 

At  midlatitudes,  magnetic  and  transverse  electric  field  mag¬ 
nitudes  at  the  surface  are  less  than  those  in  polar  regions  by  the  fac¬ 
tor  R  given  by  Equation  (35c).  Similarly,  magnitude  of  the  electro¬ 
static  field  component  at  the  surface  is  proportional  to  |B»nl,  which 
is  expressed  as  a  function  of  latitude  and  direction  of  wave  propaga¬ 
tion  by  Equation  (35b).  Polarization  of  field  profiles,  however,  also 
change  with  latitude.  We  forego  describing  polarization  changes  with 
latitude,  but  remark  that  axes  of  polarization  ellipses  are  no  longer 
aligned  vertically  and  horizontally. 


39 


IV.  ELECTROMAGNETIC  FIELDS  INDUCED  BY 
PROGRESSIVE  WAVES  IN  SHALLOW  SEAS 


In  this  section,  we  describe  electromagnetic  fields  produced 
by  progressive  waves  in  a  sea  that  is  shallower  than  a  wavelength 
(kD<  1).  We  again  suppose  that  wave  amplitudes  are  much  smaller 
than  a  wavelength.  As  in  Section  HI,  we  first  consider  surface  waves 
on  a  homogeneous  ocean  and  then  internal  waves  propagating  on  a  sharp 
thermocline  at  a  depth  d. 

Information  on  suboceanic  conductivity  profiles  is  sparse. 
Conductivity  of  the  Earth’s  crust  below  the  oceans  is  much  less  than 
that  of  sea  water.  Conductivity  increases  with  depth  within  the  Earth’s 
mantle  below  the  crust  and  becomes  much  greater  than  that  of  sea 
water.  Estimates  of  the  extent  of  the  poorly  conducting  crustal  layer 
beneath  the  oceans  range  from  a  few  tens  of  kilometers  to  a  few  hun¬ 
dreds  of  kilometers  [Coj^e^aJL,  1970J. 

To  account  for  reflections  from  suboceanic  strata,  we  choose 
a  two  layer  model  to  describe  conductivity  profiles  beneath  the  ocean 
floor.  We  suppose  that  a  poorly  conducting  layer,  extending  to  a  depth 
D+h,  rests  on  a  highly  conducting  layer  extending  indefinitely  below 
the  poorly  conducting  layer. 

A.  SURFACE  WAVES 

We  express  the  profile  of  the  velocity  potential  for  a  progres¬ 
sive  surface  wave  as 


Preceding  page  blank 
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(45a) 


o 

$  =  i^0  (cosh  kz  -  v  sinh  kz) 

where  <PQ  is  the  value  of  the  profile  at  the  surface.  As  before,  Bernoulli's 
equation  together  with  the  kinematical  relation  give  the  relation 

7?2  =  *>2/gk  (45b) 

for  the  value  of  the  polarization  coefficient  at  the  surface.  Because 
the  vertical  component  of  velocity  vanishes  at  the  ocean  floor, 

V2  =  tanhkD  ,  (45c) 

which  is  the  dispression  relation  for  surface  waves  on  an  ocean  of 
depth  D. 

We  express  the  profile  of  the  electromagnetic  potential  within 
the  sea  as  the  sum  of  exponentials  given  by  Equation  (39a).  Conse¬ 
quently,  electric  and  magnetic  fields  within  the  sea  are  expressed  by 
Equations  (40a)  and  (40b)  with  V  given  by  Equation  (45c).  As  always, 
field  profiles  above  the  surface  are  expressed  by  Equations  (25a)  and  (25b). 

The  velocity  potential  vanishes  in  suboeeanic  strata,  and  the 
magnetic  field  profile  is  expressed  by 
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in  the  poorly  conducting  layer  (D<z<D+h),  and  by 


in  the  highly  conducting  layer,  z  >D+h.  The  complex  length  h  is  deter¬ 
mined  by  the  relation 


-2o^(h-h) 

e 


~  al 
+  al 


(46c) 


where  corresponds  to  the  conductivity  of  the  poorly  conducting  layer 
and  otg,  to  the  conductivity  of  the  highly  conducting  layer.  The  electric 
field  in  suboceanic  strata  is  a  transverse  electric  type  field,  and  its 
profile  in  each  layer  is  expressed  as 


We  note  that  as  conductivity  of  the  highly  conducting  layer 
becomes  indefinitely  large  and  that  of  the  poorly  conducting  layer 
becomes  indefinitely  small,  the  length  K  approaches  h  and  approaches 
k.  Perfect  reflection  then  occurs  at  the  highly  conducting  layer,  and 
losses  within  the  poorly  conducting  layer  vanish.  Furthermore,  for 
ocean  wavelengths  that  are  much  less  than  the  extent  of  the  poorly  con¬ 
ducting  layer  (kh»l),  conductivity  profiles  of  suboceanic  strata  are 


well  approximated  by  a  vacuum  region  of  indefinite  extent  below  the 
ocean  floor. 

By  requiring  continuity  of  horizontal  components  of  electric 
and  magnetic  fields  at  the  surface  and  bottom  of  the  sea,  we  obtain 
expressions  for  the  constants  determining  electromagnetic  potentials. 
Complete  expressions  for  the  electromagnetic  potential  in  each  layer 
are  listed  in  the  Appendix.  Here,  we  express  fields  to  first  order  in 
S2. 

From  Equations  (42a)  and  (45c),  we  find  that  the  square  of  the 
ratio  of  ocean  wavelength  tc  electromagnetic  wavelength  is  expressed 
in  terms  of  wave  number  of  a  surface  wave  as 


_r-\  tanh  kD 

gD-H5- 


1/2 


(47) 


o 

The  quantity  8/2  is  shown  as  a  function  of  kD  and  DA  in  Figure  9  for 

2 

four  ocean  depths.  As  we  noted  in  Section  HI,  8  /2  is  always  small 

2 

for  wavelengths  less  than  an  ocean  depth.  We  also  observe  that  6  /2 
is  small  for  wavelengths  ten  times  an  ocean  depth  provided  ocean  depths 
are  less  than  a  kilometer,  which  we  refer  to  as  shallow  seas. 

For  shallow  seas  then,  we  again  express  the  ratio  a/k  to  first 
o 

order  in  8  as  given  by  Equation  (30d).  Furthermore,  because  the 
thickness,  h,  of  the  low  conductivity  layer  beneath  the  oceans  exceeds 
a  few  tens  of  kilometers,  we  represent  strata  below  shallow  seas  as 
a  vacuum  region  of  indefinite  extent. 

We  express  magnetic  field  profiles  in  shallow  seas  to  first 
2 

order  in  8  ,  then  as 
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k  -  -%  4  [*•*>  ♦  ^  H 


e*Zu*  (48a) 


above  the  surface  (z<0),  and  as 


■ 


•*<„  kS! 
sinh  kD  2 


(U‘0*)  e"*5^  sinh  kz  +  (H*0)  kz 


(B‘0)  sinh  k(D-z)  +  (B*0*)k(D-z)e”^^"z^ 


(48b) 


within  the  sea  (0<z<D).  Corresponding  electric  field  profiles  are  given 
by 


?a  -  -  »f0(5-fi)  ekz/2C* 


above  the  surface  and  by 
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**-0  (B-n)  r  k(D-z)  a 
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(48c) 


+  e-k(D-z) 
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within  the  sea. 

Magnitudes  of  electric  and  magnetic  fields  produced  at  the  sur¬ 
face  and  sea  bottom  by  a  one  meter  high  surface  wave  propagating  in 
polar  regions  are  shown  in  Figure  10  as  functions  of  kD  for  a  sea  depth 
of  one  kilometer.  Solid  curves  correspond  to  magnetic  field  magnitudes 
in  gamma/m  given  on  the  left-hand  coordinate  scale,  and  dashed  curves, 
to  electric  field  magnitudes  in  /*V/m  given  on  the  right-hand  coordinate. 
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For  a  fixed  value  of  kD,  field  magnitudes  are  proportional  to  / D. 

As  always,  field  magnitudes  increase  with  decreasing  wave  number. 
Eventually,  field  magnitudes  decrease  with  decreasing  wave  number 
as  ocean  wavelengths  become  comparable  to  electromagnetic  wave¬ 
lengths. 

Magnitude  and  polarization  of  magnetic  field  profiles  in  polar 
regions  are  shown  as  functions  of  z/D  in  Figure  11  for  kD  =  0.1  and  a 
sea  depth  of  one  kilometer.  We  observe  that  magnetic  field  profiles 
produced  by  ocean  waves  with  wavelengths  much  greater  than  a  sea 
depth  are  symmetrical.  Magnitudes  at  the  surface  and  sea  bed  are 
equal,  and  a  relative  minimum  occurs  at  a  depth  of  D/2.  The  field  is 
circularly  polarized  in  a  right-handed  sense  at  and  above  the  surface, 
linearly  polarized  with  a  vertical  alignment  at  a  depth  of  D/2,  and  cir¬ 
cularly  polarized  with  a  left-handed  sense  of  rotation  at  and  below  the 
sea  bed. 

Fields  induced  in  shallow  seas  in  equatorial  regions  differ 
markedly  from  those  in  polar  regions.  Field  magnitudes  in  equatorial 
regions  not  only  are  less  than  those  in  polar  regions  by  a  factor  of 
cos  6/2,  as  before,  but  also  decrease  with  wave  number  at  small  wave 
numbers,  as  shown  by  the  solid  curves  in  Figure  1 2,  which  displays 
magnetic  field  magnitudes  at  the  ocean  surface  and  bottom  that  are 
produced  by  one  meter  high  surface  waves  propagating  along  magnetic 
meridians  (6=0)  in  equatorial  regions.  Dashed  curves  in  Figure  12 
delineate  magnitudes  of  the  electrostatic  component  of  the  electric  field 
at  the  ocean  surface  and  bottom  produced  by  one  meter  high  surface 
waves  propagating  across  magnetic  meridians  (6= n/2)  in  equatorial 
regions. 
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Profiles  of  magnetic  field  magnitude  and  polarization  gen¬ 
erated  by  a  one  meter  high  surface  wave  propagating  along  magnetic 
meridians  in  equatorial  regions  are  shown  in  Figure  13a,  which  cor¬ 
responds  to  an  ocean  depth  of  1  km  and  kD  =  0.1 .  We  observe  that 
profiles  are  no  longer  symmetrical,  as  in  polar  regions,  and  that  a 
relative  minimum  occurs  at  a  depth  of  D/4  rather  than  D/2.  Elonga¬ 
tion  of  polarization  ellipses,  however,  retains  a  vertical  bias. 

Magnitude  and  polarization  profiles  of  the  electrostatic  part 
of  the  field  are  shown  as  functions  of  z/D  in  Figure  13b.  Profiles 
correspond  to  a  one  meter  high  surface  wave  propagating  across  mag¬ 
netic  meridians  in  equatorial  regions  and  having  a  wavelength  corres¬ 
ponding  to  kD  =  0.1  in  an  ocean  1  km  deep.  We  observe  that  magnitude 
of  the  field  within  the  sea  is  effectively  independent  of  depth  and  that 
its  polarization  is  highly  elongated  with  a  vertical  alignment.  At  the 
ocean  floor,  the  field  is  linearly  polarized  in  a  vertical  direction  and 
so  vanishes  in  suboceanic  strata.  Above  the  surface,  however,  the 
field  is  circularly  polarized,  and  its  magnitude  is  equal  to  the  magni¬ 
tude  of  the  horizontal  field  component  immediately  below  the  surface. 

At  intermediate  latitudes,  field  magnitudes  are  less  than 
those  in  polar  regions  by  the  factor  R  given  by  Equation  (35c),  and  the 
decrease  in  field  magnitudes  at  small  wave  numbers  lessens. 

B.  INTERNAL  WAVF 6 

For  an  ocean  of  depth  D,  we  write  the  profile  of  the  velocity 
potential  for  a  progressive  internal  wave  as 

t/>2  =  ^  [cosh  k(z-d)  -  pd  sinhk(z-d)]  (49a) 
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below  the  thermocline  (d<z<D).  The  profile  above  the  thermocline, 
is  expressed  by  Equation  (36a).  The  vertical  component  of  velocity 
vanishes  at  the  ocean  floor,  so  that 


p^  =  tanh  k(D-d) 


(49b) 


Continuity  of  pressure  and  the  vertical  component  of  velocity  at  the 
thermocline  interface  then  gives  the  dispersion  relation 


rj2  Ap/p _ 

coth  kd  +  coth  k(D-d) 


(49c) 


2  2 

where  i?  =  <u  /gk  as  before.  Continuity  of  the  vertical  component  of 
velocity  alone  at  the  thermocline  interface  gives  the  relation 


-Pd^d 


sinh  kd  (1 


V2  coth  kd) 


(49d) 


between  values  of  respective  profiles  at  the  surface,  <JfQ,  and  the  thermo¬ 
cline  interface,  i!/^.  Finally,  the  kinematical  relation  gives  the  expres¬ 
sion 


(49e) 


for  the  thermocline  displacement, 

As  before,  electric  and  magnetic  field  profiles  above  the  sur¬ 
face  are  expressed  by  Equations  (25a)  and  (25b),  and  fields  above  the 
thermocline,  by  Equations  (40a)  and  (40b).  Below  the  thermocline,  we 
express  the  profile  of  the  electromagnetic  potential  as  the  sum  of  expo¬ 
nentials 
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*s2  ~  C3  e'“<D"Z)  +  C4  ^D" Z> 


(50a) 


and  corresponding  electric  and  magnetic  field  profiles  as 


es2  =  *1**2  + 
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We  again  represent  suboceanic  strata  as  a  vacuum  region  of  indefinite 
extent  and  so  express  electric  and  magnetic  field  profiles  below  the 
ocean  as 
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(50d) 
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and 


(eK  ■  ^  2C5e-k(z-D)0 


(50e) 


Constants  specifying  the  electromagnetic  potential  profile  in 
each  layer,  obtained  by  requiring  continuity  of  horizontal  field  vectors 
at  each  interface,  are  listed  in  the  Appendix.  Here,  we  express  field 
profiles  to  first  order  in  82. 

By  using  Equation  (49c),  we  express  the  square  of  the  ratio 
of  ocean  wavelength  to  electromagnetic  wavelength  as 


so  that 


2k 


dg  (AP/P) 

kd  (coth  kd  +  coth  k(D-d) 


(1/2 


82 

T 


(4  X  l(f 7)  x  [dg  M°)(  1  -  d/D)]1/2 


(51a) 


(51b) 


In  both  deep  and  shallow  seas,  then,  internal  wavelengths  less  than  a 
few  hundred  kilometers  are  much  smaller  than  corresponding  electro¬ 
magnetic  wavelengths.* 


*  For  internal  wavelengths  of  the  order  of  100  km  or  greater,  how¬ 
ever,  internal  wave  periods  are  comparable  to  tidal  periods,  so  that 
Coriolis  forces  influence  wave  propagation  and  flow  fields  are  no 
longer  irrotational. 


50 


2 

By  neglecting  v  compared  to  unity,  as  before,  and  using  the 
first  order  expression  for  the  ratio  o/k,  as  given  by  Equation  (30d),  w 
express  field  profiles  to  first  order  in  82  as 


(52a) 


ik£d52  sinh  kD  ^  e«kd 
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above  the  surface  (z<0),  as 
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above  the  thermocline  (0<z<d),  and  as 
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(52f) 


below  the  thermocline  (d<z<D). 

For  purposes  of  illustration,  we  consider  fields  produced  by 
a  one  meter  high  internal  wave  propagating  along  a  100  m  deep  thermo¬ 
cline  in  an  ocean  1  km  deep.  We  first  consider  fields  induced  in  polar 
regions  and  then  equatorial  regions. 

Magnitudes  of  electric  and  magnetic  fields  produced  at  the 
surface,  bottom,  and  at  the  thermocline  interface  in  polar  regions 
are  shown  as  functions  of  kD  in  Figure  14.  Solid  curves  delineate 
magnetic  field  strengths  per  unit  wave  height  in  gamma/m  correspond¬ 
ing  to  the  left-hand  coordinate  scale,  and  dashed  curves,  electric  field 
strengths  per  unit  wave  height  in  microvolts/m  corresponding  to  the 
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right-hand  coordinate  scale.  For  fixed  values  of  kD  and  the  depth  ratio 
d/D,  field  magnitudes  increase  as  the  square  root  of  thermocline  depth. 
As  for  a  deep  ocean,  field  magnitudes  at  the  surface  decrease  rapidly 
for  internal  wavelengths  less  than  a  thermocline  depth.  Furthermore, 
field  magnitudes  at  the  surface  decrease  rapidly  for  wavelengths  greater 
than  an  ocean  depth.  Consequently,  internal  waves  with  wavelengths 
in  the  range,  d<  A<D,  alone  produce  an  appreciable  field  at  the  surface. 

Profiles  of  magnetic  field  magnitude  and  polarization  produced 
in  polar  regions  by  internal  waves  having  wavelengths  much  greater 
than  an  ocean  depth,  kD  =  0.5,  are  shown  in  Figure  15.  Except  near 
the  ocean  floor,  profiles  are  similar  to  those  in  a  deep  ocean,  shown 
in  Figure  4.  An  additional  polarization  reversal  occurs  near  the  ocean 
floor.  For  internal  wavelengths  much  less  than  an  ocean  depth  but  yet 
much  greater  than  a  thermocline  depth,  we  obtain  profiles  as  shown  in 
Figure  4.  Finally,  for  internal  wavelengths  much  less  than  a  thermo¬ 
cline  depth,  and  so  also  much  less  than  an  ocean  depth,  we  obtain  pro¬ 
files  as  shown  in  Figure  5 . 

In  equatorial  regions,  magnitudes  of  the  transverse  electric 
part  of  the  electric  field  and  the  magnetic  field  are  again  less  than 
values  in  polar  regions  by  a  factor  of  cos  0/2.  Furthermore,  polari¬ 
zation  profiles  and  the  dependence  of  field  magnitudes  at  the  thermo¬ 
cline  interface  on  wave  number  change.  Magnitudes  of  the  magnetic 
field  and  the  electrostatic  part  of  the  electric  field  in  equatorial  regions 
are  shown  as  functions  of  kD  in  Figure  16. 

Solid  curves  delineate  magnetic  field  strengths  produced  at 
the  surface,  bottom,  and  at  the  thermocline  interface  per  unit  wave 
height  in  gamma/m,  corresponding  to  the  left-hand  coordinate  scale, 
by  internal  waves  propagating  along  magnetic  meridians  (0=0).  Dashed 


curves  represent  magnitudes  of  the  electrostatic  part  of  the  electric 

field  produced  per  unit  wave  height  immediately  below  the  surface, 

f®o+l/^d’  immediately  above  and  below  the  thermoeline  interface, 

and  \E^+\/(^  respectively,  and  immediately  above  the  ocean 

floor  Electric  field  magnitudes  correspond  to  internal  waves 

propagating  across  magnetic  meridians  {0  =  n/2)  and  are  expressed  in 
o 

microvolts/m  ,  corresponding  to  the  right-hand  coordinate  scale. 

We  note  that  the  magnetic  field  strength  at  the  thermoeline  interface 
now  also  decreases  for  small  wave  numbers. 

Profiles  of  magnetic  field  magnitude  and  polarization  in  equa¬ 
torial  regions  are  shown  in  Figure  17a  for  kD  =  0. 5  and  0=0.  We 
observe  that  polarization  ellipses  are  predominantly  aligned  vertically 
as  opposed  to  a  predominant  horizontal  alignment  in  polar  regions. 

Finally,  profiles  of  magnitude  and  polarization  of  the  electro¬ 
static  part  of  the  electric  field  are  shown  in  Figure  17b.  We  observe 
that  the  field  is  linearly  polarized  with  a  vertical  alignment  at  both 
the  surface  and  ocean  floor  and  that  the  electrostatic  field  vanishes 
above  and  below  the  sea. 
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V.  CONCLUSION 


To  conclude,  we  summarize  salient  features  of  the  electro¬ 
magnetic  field  induced  by  progressive  ocean  waves  and  mark  features 
germane  to  spectral  measurements.  First,  the  electromagnetic  field 
is  the  sum  of  a  transverse-electric  type  field  and  an  electrostatic  type 
field. 

The  electrostatic  part  of  the  field  is  polarized  in  a  vertical 
plane  containing  the  direction  of  wave  propagation.  Polarization  of 
the  field  is  independent  of  latitude.  At  the  ocean  floor,  polarization 
is  linear  with  a  vertical  alignment,  for  both  surface  and  internal  waves, 
so  that  the  electrostatic  field  vanishes  in  suboceanic  strata.  Within 
the  sea,  polarization  of  the  electrostatic  field  produced  by  surface 
waves  is  elliptical  with  a  vertical  major  axis  and  a  left-handed  sense 
of  rotation.  Polarization  ellipses  are  highly  elongated  for  wavelengths 
of  surface  waves  that  are  much  greater  than  an  ocean  depth  and  are 
nearly  circular  for  wavelengths  much  less  than  ocean  depth.  Above 
the  surface,  polarization  is  circular  with  a  right-handed  sense  of  rota¬ 
tion. 

Polarization  of  the  electrostatic  'ield  produced  by  internal 
waves,  however,  is  linear  with  a  vertical  alignment  at  the  surface,  so 
that  the  electrostatic  part  of  the  field  vanishes  above  the  surface.* 


*  A  vertically  aligned  linear  polarization  at  the  surface  results  because 
surface  displacements  produced  by  internal  waves  are  negligibly  small. 
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Within  the  sea,  polarization  is  elliptical  with  a  vertical  major  axis 
and  a  right-handed  sense  of  rotation  above  the  thermocline  and  a  left- 
handed  sense  of  rotation  below  the  thermocline.  Polarization  ellipses 
above  the  thermocline  are  highly  elongated  for  wavelengths  of  internal 
waves  that  are  much  greater  than  a  thermocline  depth  and  are  nearly 
circular  (except  near  the  surface)  for  wavelengths  much  greater  than 
a  thermocline  depth.  Below  the  thermocline,  polarization  ellipses 
are  highly  elongated  for  wavelengths  much  less  than  an  ocean  depth 
and  are  nearly  circular  (except  near  the  ocean  floor)  for  wavelengths 
much  less  than  an  ocean  depth. 

Magnitude  of  the  electrostatic  part  of  the  field  is  proportional 
to  the  horizontal  component  of  the  Earth's  magnetic  field  that  is  trans¬ 
verse  to  the  direction  of  wave  propagation  and  so  vanishes  in  polar 
regions  and  for  waves  propagating  along  magnetic  meridians.  Magni¬ 
tudes  are  greatest  in  equatorial  regions  for  waves  propagating  across 
magnetic  meridians.  Maximum  magnitudes  of  the  electrostatic  field 
immediately  below  the  surface  produced  per  unit  surface  wave  height 
are  a  few  microvolts/m  for  wavelengths  much  greater  than  an  ocean 
depth  and  increase  as  the  square  root  of  wave  number  for  wavelengths 
much  less  than  an  ocean  depth.  Electrostatic  field  strength  and  cur¬ 
rent  density  at  the  surface  are  the  only  surface  quantities  that  increase 
with  wave  number. 

Maximum  electrostatic  field  strengths  produced  immediately 
below  the  surface  per  unit  internal  wave  height  are  a  few  tenths  of 
microvolts/m  for  wavelengths  much  greater  than  a  thermocline  depth 
and  decrease  rapidly  for  wavelengths  much  less  than  a  thermocline 
depth.  Electrostatic  field  strengths  at  the  thermocline  interface,  how¬ 
ever,  increase  as  the  square  root  of  wave  number. 
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The  magnetic  field  of  the  transverse  electric  part  of  the  field 
is  polarized  in  a  vertical  plane  containing  the  direction  of  wave  propa¬ 
gation,  and  the  electric  field,  by  definition,  is  linearly  polarized  in 
the  horizontal  direction  transverse  to  the  direction  of  wave  propaga¬ 
tion.  Within  the  sea,  polarization  of  the  magnetic  field  produced  by 
both  surface  and  internal  waves  varies  markedly  with  both  depth  and 
latitude.  At  and  above  the  surface,  however,  polarization  is  always 
circular  with  a  right-handed  sense  of  rotation,  whatever  the  type  of 
wave  motion  or  latitude. 

Magnitude  of  the  magnetic  field  is  proportional  to  the  product 
of  wave  height'  and  phase  speed  of  ocean  waves,  provided  wavelengths 
are  smaller  than  electromagnetic  wavelengths  and  much  greater  than 
wave  heights.  Magnitude  of  the  transverse  electric  field  is  the  product 
of  phase  speed  and  the  vertical  component  of  the  induced  magnetic  field 
and  thus  is  proportional  to  the  product  of  wave  height  and  the  square 
of  phase  speed.  As  a  result,  field  magnitudes  per  unit  wave  height 
increase  with  decreasing  wave  number  and  attain  maximum  values  at 
small  wave  numbers  —  much  as  do  phase  speeds.  Magnitudes  are 
also  greater  in  polar  regions  than  equatorial  regions  because  the  Earth's 
field  is  stronger  in  polar  regions. 

In  polar  regions,  magnetic  and  electric  field  strengths  at  the 
ocean  surface  and  bottom  produced  by  surface  waves  are  reinforced 
by  surface  and  bottom  reflections  for  wavelengths  greater  than  an  ocean 
depth.  Similarly,  field  strengths  at  a  thermocline  interface  produced 
by  internal  waves  are  reinforced  by  reflections,  but  field  strengths  at 
the  surface  and  bottom  are  reduced  by  interfering  reflections  at  wavelengths 


*  Thermocline  displacement  is  the  wave  height  of  an  internal  wave. 
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greater  than  an  ocean  depth.  Surface  waves  produce  maximum  mag¬ 
netic  and  electric  field  magnitudes  per  unit  wave  height  of  the  order 
of  10  gamma/m  and  1  microvolt/m  ,  respectively,  at  the  surface  in 
polar  regions.  Maximum  field  magnitudes  per  unit  wave  height  pro¬ 
duced  at  the  surface  by  internal  waves  are  of  the  order  of  0.1  gamma/ 
m  and  10"^  microvolts/m^  in  polar  regions. 

In  equatorial  regions,  however,  surface  and  bottom  reflec¬ 
tions  interfere  and  so  reduce  magnetic  and  electric  field  strengths  at 
the  ocean  surface  and  bottom  produced  by  surface  waves  with  wave¬ 
lengths  much  greater  than  an  ocean  depth,  as  well  as  field  strengths 
at  a  thermocline  interface  produced  by  internal  waves.  Maximum  field 
strengths  in  equatorial  regions  are  about  one  half  those  in  polar  regions, 
for  waves  propagating  along  magnetic  meridians.  The  transverse  elec¬ 
tric  part  of  the  field  vanishes  for  waves  propagating  across  magnetic 
meridians  in  equatorial  regions. 

For  purposes  of  spectral  measurements,  measurement  of  the 
induced  magnetic  field  at  or  above  the  surface  is  most  appropriate. 
Because  electric  field  strengths  are  small,  electrode  separations  of 
the  order  of  tens  of  meters  are  required,  which  complicates  interpre¬ 
tation  of  spectral  measurements.  Magnetic  field  measurements,  how¬ 
ever,  are  effectively  point  measurements,  so  that  spectral  features  of 
the  induced  magnetic  field  are  simply  related  to  those  of  ocean  waves. 
The  power  spectrum  of  the  induced  magnetic  field  is  simply  the  pro¬ 
duct  of  the  ocean  wave  power  spectrum  and  the  squared  modulus  of  a 
transfer  function,  which  is  the  magnetic  field  profile  induced  per  unit 
wave  height  by  a  progressive  wave.  Consequently,  dependence  of  the 
magnetic  field  profile  induced  per  unit  wave  height  at  or  above  the  sur¬ 
face  on  wave  number  is  the  feature  of  the  electromagnetic  field  induced 
by  a  progressive  ocean  wave  that  is  germane  to  spectral  measurement. 
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In  this  regard,  we  note  that  magnetic  field  spectra  are  softer 
than  ocean  wave  spectra  because  induced  fields  are  greater  at  small 
wave  numbers.  Furthermore,  spectra  soften  with  increasing  altitude, 
as  the  exponential  decrease  in  field  magnitudes  with  altitude  is  lesser 
at  smaller  wave  numbers .  We  defer  detailed  description  of  spectral 
features,  however,  to  a  subsequent  report. 
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FIGURE  CAPTIONS 


Figure  1 .  Magnitudes  of  electric  and  magnetic  fields  produced  at  the 
surface  of  a  deep  ocean  by  a  one  meter  high  surface  wave  propagating 
in  polar  regions  shown  as  functions  of  wave  number  k  along  the  lower 
abscissa  and  as  functions  of  reciprocal  wavelength,  lA,  along  the 
upper  abscissa.  The  solid  curve  delineates  magnetic  field  strength 
per  unit  wave  height  in  gamma/m  and  corresponds  to  the  left-hand 
coordinate  scale.  The  dashed  curve  delineates  electric  field  strength 
per  unit  wave  height  in  microvolts/m  and  corresponds  to  the  right- 
hand  coordinate  scale.  Straight  dashed  lines  are  asymptotes  for  large 
and  small  wave  numbers. 

Figure  2.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  a  deep  ocean  by  a  one  meter  high  surface  wave  propagating 
in  polar  regions  with  a  wavelength  of  10  m  corresponding  to  a  frequency 
of  about  0.4  Hz.  Altitude  and  depth  are  measured  in  units  of  an  ocean 
wavelength. 

Figure  3.  Magnitudes  of  electric  and  magnetic  fields  produced  at  the 
surface  and  at  a  100  meter  deep  thermocline  by  a  one  meter  high  inter¬ 
nal  wave  propagating  in  a  deep  ocean  in  polar  regions  shown  as  func¬ 
tions  of  dimensionless  wave  number  kd  along  the  lower  abscissa  and 
as  functions  of  dA  along  the  upper  abscissa,  where  d  denotes  thermo¬ 
cline  depth.  Solid  curves  delineate  magnetic  field  strengths  per  unit 
wave  height  in  gamma/m  and  correspond  to  the  left-hand  coordinate 
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scale.  Dashed  cu  *ves  delineate  electric  field  strengths  per  unit  wave 
height  in  mi^rovolts/m  and  correspond  to  the  right-hand  coordinate 
scale. 

Figure  4.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  a  deep  ocean  containing  a  sharp  thermocline  at  a  depth  of 
100  meters  by  a  one  meter  high  internal  wave  propagating  in  polar 
regions  and  having  a  wavelength  much  greater  than  a  thermocline  depth, 
kd  =  0. 1 .  Altitude  and  depth  are  measured  in  units  of  a  thermocline 
depth,  d. 

Figure  5.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  a  deep  ocean  containing  a  sharp  thermocline  at  a  depth  of  100 
meters  by  a  one  meter  high  internal  wave  propagating  in  polar  regions 
and  having  a  wavelength  much  less  than  a  thermocline  depth,  kd  =  4. 
Altitude  and  depth  are  measured  in  units  of  a  thermocline  depth,  d. 

Figure  6 .  Magnitudes  of  the  magnetic  field  and  the  electrostatic  part 
of  the  electric  field  produced  at  the  surface  and  at  a  100  meter  deep 
thermocline  by  a  one  meter  high  internal  wave  propagating  in  a  deep 
ocean  in  equatorial  regions  shown  as  functions  of  kd  along  the  lower 
abscissa  and  as  functions  of  d/A  along  the  upper  abscissa.  Magnetic 
field  magnitudes  correspond  to  internal  wave  propagation  along  mag¬ 
netic  meridians,  and  electric  field  magnitudes,  to  propagation  across 
magnetic  meridians.  Solid  curves  delineate  magnetic  field  strengths 
per  unit  wave  height  in  gamma/m  and  correspond  to  the  left-hand 
coordinate  scale .  Dashed  curves  delineate  electric  field  strengths 
per  unit  wave  height  in  microvolts/m  and  correspond  to  the  right- 
hand  coordinate  scale. 
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Figure  7a.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  a  deep  ocean  containing  a  sharp  thermocline  at  a  depth  of 
100  meters  by  a  one  meter  high  internal  wave  propagating  along  mag¬ 
netic  meridians  in  equatorial  regions  and  having  a  wavelength  much 
greater  than  a  thermocline  depth,  kd  =  0.1 .  Altitude  and  depth  are 
measured  in  units  of  a  thermocline  depth,  d. 

Figure  7b.  Profile  of  magnitude  and  polarization  of  the  electrostatic 
part  of  the  electric  field  produced  in  a  deep  ocean  containing  a  sharp 
thermocline  at  a  depth  of  100  meters  by  a  one  meter  high  internal 
wave  propagating  across  magnetic  meridians  in  equatorial  regions  and 
having  a  wavelength  much  greater  than  a  thermocline  depth,  kd  =  0.1. 
Altitude  and  depth  are  measured  in  units  of  a  thermocline  depth,  d. 

Figure  8a.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  a  deep  ocean  containing  a  sharp  thermocline  at  a  depth  of 
100  meters  by  a  one  meter  high  internal  wave  propagating  along  mag¬ 
netic  meridians  in  equatorial  regions  and  having  a  wavelength  much 
less  than  a  thermocline  depth,  kd  =  4.  Altitude  and  depth  are  mea¬ 
sured  in  units  of  a  thermocline  depth,  d. 

Figure  8b.  Profile  of  magnitude  and  polarization  of  the  electrostatic 
part  of  the  electric  field  produced  in  a  deep  ocean  containing  a  sharp 
thermocline  at  a  depth  of  100  meters  by  a  one  meter  high  internal 
wave  propagating  across  magnetic  meridians  in  equatorial  regions 
and  having  a  wavelength  much  less  than  a  thermocline,  kd  =  4.  Altitude 
and  depth  are  measured  in  units  of  a  thermocline  depth,  d. 
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Figure  f) .  Square  of  the  ratio  of  wavelength  of  a  surface  wave  on  an 
ocean  oi  depth  D  to  wavelength  of  an  electromagnetic  wave,  8  /2, 
shown  as  a  function  of  kD  on  the  lower  abscissa  and  as  a  function  of 
D/x  on  the  upper  abscissa  for  ocean  depths  of  10,  1,  0.5,  and  0.1 
kilometers.  Straight  dashed  lines  are  asymptotes  for  large  and 
small  values  of  kD. 

Figure  10.  Magnitudes  of  electric  and  magnetic  fields  produced  at  the 
surface  and  bottom  of  an  ocean  1  km  in  depth  by  a  one  meter  high  sur¬ 
face  wave  propagating  in  polar  regions  shown  as  functions  of  kD  along 
the  lower  abscissa  and  as  functions  of  D/x  along  the  upper  abscissa, 
where  D  denotes  ocean  depth.  Solid  curves  delineate  magnetic  field 
strengths  per  unit  wave  height  in  gamma/m  and  correspond  to  the  left- 
hand  coordinate  scale.  Dashed  curves  delineate  electric  field  strengths 
per  unit  wave  height  in  microvolts/m  and  correspond  to  the  right-hand 
coordinate  scale. 

Figure  11 .  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  an  ocean  1  km  deep  by  a  one  meter  high  surface  wave  propa¬ 
gating  in  polar  regions  and  having  a  wavelength  much  greater  than  an 
ocean  depth,  kD  =  C.l .  Altitude  and  depth  are  measured  in  units  of 
an  ocean  depth,  D. 

Figure  12.  Magnitudes  of  the  magnetic  field  and  the  electrostatic  part 
of  the  electric  field  produced  ai  the  surface  and  bottom  of  an  ocean 
1  km  in  depth  by  a  one  meter  high  surface  wave  propagating  in  equa¬ 
torial  regions  shown  as  functions  of  kD  along  the  lower  abscissa  and 
as  functions  of  D/\  along  the  upper  abscissa.  Magnetic  field  magni¬ 
tudes  correspond  to  surface  wave  propagation  along  magnetic  meridians, 
and  electric  field  magiitudes,  to  propagation  across  magnetic  meridians. 


Solid  curves  delineate  magnetic  field  strengths  per  unit  wave  height 
in  gamma/m  and  correspond  to  the  left-hand  coordinate  scale.  Dashed 
curves  delineate  electric  field  strengths  per  unit  wave  height  in  micro- 
volts/m  and  correspond  to  the  right-hand  coordinate  scale. 

Figure  13a.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  an  ocean  1  km  deep  by  a  one  meter  high  surface  wave  propa¬ 
gating  along  magnetic  meridians  in  equatorial  regions  and  having  a 
wavelength  much  greater  than  an  ocean  depth,  kD  =  0.1 .  Altitude  and 
depth  are  measured  in  units  of  an  ocean  depth,  D. 

Figure  13b.  Profile  of  magnitude  and  polarization  of  the  electrostatic 
part  of  the  electric  field  produced  in  an  ocean  1  km  deep  by  a  one  meter 
high  surface  wave  propagating  across  magnetic  meridians  in  equatorial 
regions  and  having  a  wavelength  much  greater  than  an  ocean  depth, 
kD  =  0.1.  Altitude  and  depth  are  measured  in  units  of  an  ocean  depth, 
D. 

Figure  14.  Magnitudes  of  electric  and  magnetic  fields  produced  at  the 
surface,  bottom,  and  at  a  100  meter  deep  thermocline  in  an  ocean 
1  km  deep  by  a  one  meter  high  internal  wave  propagating  in  polar 
regions  shown  as  functions  of  kD  along  the  lower  abscissa  and  as 
functions  of  D/x  along  the  upper  abscissa,  where  D  denotes  ocean 
depth.  Solid  curves  delineate  magnetic  field  strengths  per  unit  wave 
height  in  gamma/m  and  correspond  to  the  left-hand  coordinate  scale. 
Dashed  curves  delineate  electric  field  strengths  per  unit  wave  height 
in  microvolts/m  and  correspond  to  the  right-hand  coordinate  scale. 
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Figure  15.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  an  ocean  1  km  deep  containing  a  sharp  thermocline  at  a  depth 
of  100  m  by  a  one  meter  high  internal  wave  propagating  in  polar  regions 
and  having  a  wavelength  much  greater  than  an  ocean  depth,  kD  =  0.5. 
Altitude  and  depth  are  measured  in  units  of  a  thermocline  depth,  d. 

Figure  16.  Magnitudes  of  the  magnetic  field  and  electrostatic  part  of 
the  electric  field  produced  at  the  surface,  bottom,  and  at  a  100  meter 
deep  thermocline  in  an  ocean  1  km  deep  by  a  one  meter  high  internal 
wave  propagating  in  equatorial  regions  shown  as  functions  of  kD  along 
the  lower  abscissa  and  as  functions  of  D/X  along  the  upper  abscissa. 
Magnetic  field  magnitudes  correspond  to  internal  wave  propagation 
along  magnetic  meridians,  and  electric  field  magnitudes,  to  propaga¬ 
tion  across  magnetic  meridians.  Solid  curves  delineate  magnetic  field 
strengths  per  unit  wave  height  in  gamma/m  and  correspond  to  the  left- 
hand  coordinate  scale.  Dashed  curves  delineate  electric  field  strengths 
per  unit  wave  height  in  microvolts/m  and  correspond  to  the  right-hand 
coordinate  scale. 

Figure  17a.  Profile  of  magnetic  field  magnitude  and  polarization  pro¬ 
duced  in  an  ocean  1  km  deep  containing  a  sharp  thermocline  at  a  depth 
of  100  m  by  a  one  meter  high  internal  wave  propagating  along  magnetic 
meridians  in  equatorial  regions  and  having  a  wavelength  much  greater 
than  an  ocean  depth,  kD  =  0.5,  Altitude  and  depth  are  measured  in 
units  of  a  thermocline  depth,  d. 

Figure  17b.  Profile  of  magnitude  and  polarization  of  the  electrostatic 
part  of  the  electric  field  produced  in  an  ocean  1  km  deep  containing  a 
sharp  thermocline  at  a  depth  of  100  m  by  a  one  meter  high  internal 
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wave  propagating  across  magnetic  meridians  in  equatorial  regions 
and  having  a  wavelength  much  greater  than  an  ocean  depth,  kD  =  0.5. 
Altitude  and  depth  are  measured  in  units  of  a  thermocline  depth,  d. 
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APPENDIX.  PROFILES  OF  ELECTROMAGNETIC  POTENTIALS 


Profiles  of  electromagnetic  fields  induced  by  progressive  waves 

in  horizontally  stratified  oceans  are  specified  by  an  electromagnetic 

potential  profile,  ,  and  a  velocity  potential  profile,  ^ ,  as  expressed 

s 

by  Equations  (20a)  and  (20b)  in  Section  EL  The  transverse-electric 
part  of  the  field  depends  on  both  electromagnetic  and  velocity  poten¬ 
tial  profiles,  but  the  electrostatic  part  of  the  field  depends  on  the  velocity 
potential  profile  alone.  In  Sections  HI  and  IV,  we  described  field  pro¬ 
files  produced  in  deep  and  shallow  seas,  respectively,  by  progressive 
ocean  waves  having  wavelengths  much  less  than  wavelengths  of  cor¬ 
responding  electromagnetic  waves.  First  order  expressions  for  field 
profiles  are  quite  adequate  provided  wavelengths  of  surface  waves  are 
less  than  about  10  km  and  wavelengths  of  internal  waves  are  less  than 
about  100  km,  and  so  are  satisfactory  for  most  purposes. 

For  completeness,  however,  we  give  here  exact  expressions 
for  electromagnetic  potential  profiles  generated  by  surfaces  waves  in 
an  ocean  of  depth  D  and  by  internal  waves  in  an  ocean  of  depth  D  con¬ 
taining  a  sharp  thermocline  at  a  depth  d.  As  in  Section  IV,  we  sup^se 
that  suboceanic  conductivity  profiles  are  adequately  described  by  a 
poorly  conducting  layer,  extending  to  a  depth  D+h,  resting  on  a  highly 
conducting  layer  extending  indefinitely  below  the  poorly  conducting 
layer.  But  we  do  not  further  suppose,  as  in  Section  IV,  that  the  poorly 
conducting  layer  is  effectively  a  vacuum  region  of  indefinite  extent, 
although  the  supposition  is  adequate  for  ocean  wavelengths  less  than 
a  few  tens  of  kilometers. 
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A.  ELECTROMAGNETIC  POTENTIAL  PROFILE  GENERATED 

BY  SURFACE  WAVES 

The  profile  of  the  electromagnetic  potential  generated  by  a 
progressive  surface  wave  in  an  ocean  of  depth  D  is  given  by 

*a  *  A0ekz  (Ala) 

above  the  surface  (z<0),  by 

<t>s  =  C i  eaZ  +  C2  e"aZ  (Alb) 

within  the  sea  (0<z<D),  by 


*cl  =  2Cg  cosh  [z  -  (D+h)] 


(Ale) 


within  the  poorly  conducting  layer  (D<z<D+h),  and  by 


-«2[z-(D+h)] 

e 


within  the  highly  conducting  layer  (z>D-th),  where 


-2a1(h-h) 

e 


(Aid) 


(Ale) 


defines  the  complex  length  h  and  aj  and  correspond  to  conductivities 
of  the  poorly  conducting  and  highly  conducting  layers,  respectively. 

The  constants  Aq,  C^,  C 2,  and  C^  are  expressed  in  terms  of 
the  value  of  the  velocity  potential  profile  at  the  surface,  ,  by  the 
relations 
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V2  (f  kD  _.«l  ^1  .  . 

i'y+  cosh  kD([y-  e  "  2^T  **  1311,1  «ih)  (B-0) 

+  y+e  kD  “  +  tanh  (B-  r*)j 


Tr+cosh  kD  {[2(tt  -  tenh  “jh)  e'(a_k)D 

'  I1  +  |)(t  "  tanh  “lbj  (B'C) 

"f1  +tanh  <*!*)  (B-(N-)j 

°2  =  Ty+cosh  kD | f  2("cT_  +tanha1hje(tt+k)D 

+  f  "  Mt  "  tanh  ^b)  (5-0) 


+  181111  “ibj  <§•{>*) 


(A  2a) 


(A2b) 


(A  2c) 


2C3sinha1h  =  Ct  eaD  +  C0  e'0^  +  ^2^ 
J  1  1  2 e  cosh  kD 


where 


(A  2d) 
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B.  ELECTROMAGNETIC  POTENTIAL  PROFILE  GENERATED 

BY  INTERNAL  WAVES 

The  profile  of  the  electromagnetic  potential  generated  by  a 
progressive  internal  wave  in  an  ocean  of  depth  D  containing  a  sha^p 
thermocline  at  a  depth  d  is  given  by 

4>a  =  A()  ekz  (A4a) 

above  the  surface  (z<0),  by 

=  cn  e“z  +  ci2e‘°z  <A4b) 
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above  the  thermocline  (d<z<D),  by 

*S2  ■  e“<Z'd)  ♦  C22  s'“(Z'd)  <**> 

below  the  thermocline  (d<z<D),  by 

4>cl  =  2Cg  cosh  aj  [z  -  (D+h)]  (A4d) 


within  the  poorly  conducting  layer  (d  <  z  <  D+h),  and  by 


ai(h-h)  -a^z-CD+h)] 
e  e 


(A4e) 


within  the  h.ghly  conducting  layer  (z>D+h),  where  the  complex  length 
h  is  defined  as  before  [Equation  (Ale)]. 

The  constants  Aq,  C^,  Cgj,  C22>  and  Cg  are  expressed 
in  terms  of  the  value  of  the  velocity  potential  profile  at  the  surface,  ^  , 
by  the  relations 


2  y+  sinh  k(D-d) 


7  sinh  k(D-d)  -  G  sinh  kD 


+ 


7  sinh  k(D-d)  -  G  sinh  kD 

“r  *f* 


/« l 

+  2  sinh  kd  I  -^-  +  tanh 


(B*t)*) 


(A5a) 
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r  t 


4  |  cn  =  -V5j(f  -  l)<B-C>)  +  (a  + 1)  (B-  0*) 


4  «  r 
4  k  l12 


4  —  P 
k  L22 


and 


+  2A 


+ 1)  (B-D)  +  (f  -  lj  (B-0*) 


+  2A 


•M 


0Q  sinh  kD 
^^^21  ~  sinh  k(D-d) 


|*-l)(B-{>)  +  /g  +  lj(g>) 


+  4  a  cn  e' 


.ad 


ib  sinh  kD 
o _ 

sinh  k(D-d) 


(g+l)  (B-O)  +  (g-  l)  (B-O- 


+  4«p  e“ad 
+  *k  C12e 


(A5b) 

(A5c) 

A 

(A5d) 

)y/2 

(A5e) 
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_  -  ^  sinh  kd 

2C3sinh“ih  -arog 


+  C21  6 


a(D-d) 


(A5f) 


where 


«1  a- 

=  2  cosh  a(D-d)  y  1  *5*  tanh  a  (D-d) 
+  {%  tanh  a(D-d)  +  ij  tanh  q^K 


and  y+  is  defined  by  Equation  (A2e). 


(A5g) 


The  thermocline  displacement,  or  wave  height  of  an  internal 
wave,  q,  is  related  to  the  value  of  the  velocity  potential  profile  at  the 
surface  by  the  expression 


i  £ 


-  sinh  kd 


As  kD  becomes  large,  the  constants  C?1  and  C„  vanish, 


A  a  Jl 
0 


f  (*  "  e  (  k)d]  (hi)^  +  $•&*)  »  (A7a) 


:ilS  (*  +s)  e'(“‘k)d  (2-l|)(B.D)  +  (B-0*)  ,  (A7b) 
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(B-O*) 


(A7c) 


and 


C 


22 


,  ,  ,  2 
a  -  k\ 

a  +  k  / 


-  e"(““k>d)  (B-0*)  J  .  (A7d) 


Furthermore,  as  kd  becomes  large,  which  implies  that  kD  is 
also  large,  the  constants  C ^  Cg,  and  vanish, 


A.» 


(A8a) 


cn= 


(A8b) 


and 
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lb  e 

^ o 


kd 


2  y/2 


('  +a-)|<5^  +  (£tf)  (5-')*) 


(A8c) 
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